Hello Statalist,
I've run into an issue trying to make high-resolution (quarter-degree lat/lon) maps of temperature trends over the period from 1979 to present in the conterminous United States. I have values for each quarter-degree grid cell already calculated.
Specifically, the trends I'm working with range from 0.006 to -0.006 degrees C per month. There is considerable spatial structure to the data, so I want to use a color palette that provides me the ability to see these features. Both Reds2 and Blues2 work quite well, giving me 99 different colors to work with. BuYlRd, which combines the two, only allows 11 different colors, and produces a map that is too coarse to be of much use.
Unfortunately there are both cooling and warming trends present on the map, and using only red or only blue doesn't effectively differentiate between the two. What I want to do is create a map that combines a map of all the grid cells with warming trends (using the Reds2 color scheme) and all of the grid cells with cooling trends (using the Blues2 color scheme).
Here is the warming trend map: http://i81.photobucket.com/albums/j2...psf210b7d4.png
And the cooling trend map: http://i81.photobucket.com/albums/j2...ps65f19200.png
And the relevant code:
spmap prism using quarter_degree_cells.dta, ///
id(id) fcolor(Reds2) clnumber(98) ///
legenda(off) clmethod(custom) clbreaks(0(0.000075)0.006) ///
ocolor(none ..)
spmap prism using quarter_degree_cells.dta, ///
id(id) fcolor(Blues2) clnumber(98) ///
legenda(off) clmethod(custom) clbreaks(-0.006(0.000075)0) ///
ocolor(none ..)
Is this possible? Or am I stuck using the BuYlRd or Rainbow scheme?
Thanks,
-Zeke
I've run into an issue trying to make high-resolution (quarter-degree lat/lon) maps of temperature trends over the period from 1979 to present in the conterminous United States. I have values for each quarter-degree grid cell already calculated.
Specifically, the trends I'm working with range from 0.006 to -0.006 degrees C per month. There is considerable spatial structure to the data, so I want to use a color palette that provides me the ability to see these features. Both Reds2 and Blues2 work quite well, giving me 99 different colors to work with. BuYlRd, which combines the two, only allows 11 different colors, and produces a map that is too coarse to be of much use.
Unfortunately there are both cooling and warming trends present on the map, and using only red or only blue doesn't effectively differentiate between the two. What I want to do is create a map that combines a map of all the grid cells with warming trends (using the Reds2 color scheme) and all of the grid cells with cooling trends (using the Blues2 color scheme).
Here is the warming trend map: http://i81.photobucket.com/albums/j2...psf210b7d4.png
And the cooling trend map: http://i81.photobucket.com/albums/j2...ps65f19200.png
And the relevant code:
spmap prism using quarter_degree_cells.dta, ///
id(id) fcolor(Reds2) clnumber(98) ///
legenda(off) clmethod(custom) clbreaks(0(0.000075)0.006) ///
ocolor(none ..)
spmap prism using quarter_degree_cells.dta, ///
id(id) fcolor(Blues2) clnumber(98) ///
legenda(off) clmethod(custom) clbreaks(-0.006(0.000075)0) ///
ocolor(none ..)
Is this possible? Or am I stuck using the BuYlRd or Rainbow scheme?
Thanks,
-Zeke
