Re: st: question about spmap

[Sergiy Radyakin <serjradyakin@gmail.com>]

Hi, Ana,

Your data points may have coordinates measured differently from your map.
This happens often when the map projections are different. For example if your
map was in UTM you may use the following manual to convert it from the
northing-easting to the lat-long notation:

http://earth-info.nga.mil/GandG/publications/tm8358.2/TM8358_2.pdf

To check that, see if either the map or the point coordinates have
values measured
in 100thousands or millions (or consult the data provider). If so,
refer to the manual
above to convert your data. I prefer working in lat-long, but that's a matter of
convenience and context.

I do have a working code to convert from UTM to lat-long, so if this
is your situation,
when Kit is back, I could post it to the SSC.

Best, Sergiy Radyakin




On Wed, Aug 25, 2010 at 2:11 PM, Maria Ana Vitorino
<vitorino@wharton.upenn.edu> wrote:
> Dear StataList users,
> I have a question regarding spmap.
> If I do
> spmap using pt-c, id(id)
> then I get a simple map with all the boundaries for the different regions,
> which is ok.
> But, as soon as I add some point data (from another file) the boundaries
> disappear.
> For example, when I do:
> spmap using pt-c, id(id) point(data("points.dta") xcoord(xcoord)
> ycoord(ycoord) fcolor(blue) )
> only the points show up, not the boundaries for the regions.
>
> How can I have both the boundaries of the regions and the point data?
>
> Thanks for the help,
>
> Ana
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
>

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/