Help with combining 2 shapefiles

[Juan Esteban Berrio]

Hi, I have previosly worked with the "spmap" command, so Im familiar with it. But now im trying to do something that is a more advanced, Im trying to combine 2 shapefiles and I dont know if that is possible.

The Colombian shapefile that I have from a goverment source doesnt show the San Andrés y Providencia department (is like a state) properly so when I map something like lets say the unemployment rate it doesnt show correctly due to how small it is shown in the map .

So there are 2 options here, either there is a way to scale the polygon of San Andrés y Providencia and make it a little bigger to properly show what im trying to map, or there is way to combine 2 shapefiles (the one above but without San Andrés y Providencia) and one of San Andrés that I found from the same government source.


The only problem with this shapefile is that I dont know why the 2 islands are so far away from each other when I map the shapefile.

I have tried to add the coordinates database from the shapefile of San Andrés to the shapefile of Colombia (replacing the ones from the tiny San Andrés) but I always get a black map when I try to map the resulting database.

So if anyone could help me or point me in the right way, I'll be very greatful.

[Scott Merryman]

Here is one way, resizing and moving two of the larger islands. Using data from http://www.diva-gis.org/gdata

shp2dta using col_adm2, database(db) coordinates(coord) genid(id) replace   gencentroids(c)
use db,clear
spmap using coord, id(id) name(orig_size, replace)

sum x_c if id == 853
local xc1 = r(mean)
sum y_c if id == 853
local yc1 = r(mean)

sum x_c if id == 854
local xc2 = r(mean)
sum y_c if id == 854
local yc2 = r(mean)

use coord,clear

replace _X = 10*(_X - `xc1') + `xc1' + 3 if _ID == 853
replace _Y = 10*(_Y - `yc1') + `yc1' - 2  if _ID == 853
replace _X = 10*(_X - `xc2') + `xc2' + 3 if _ID == 854
replace _Y = 10*(_Y - `yc2') + `yc2' - 3 if _ID == 854
save,replace

use db
spmap using coord, id(id)

Original map:

Resized:

[Juan Esteban Berrio]

Scott that is just what I needed thank you very much for your help.

I would try that with my shapefile.