I have a data set in one projection, and a shapefile in another, and I am trying to combine the two. I am using Python. I apologize for posting a "help me transform coordinates" question on this site - there are so many.
I assume the way to do this conversion is with proj4. I am unable to find/understand the proj4 string needed to do this conversion.
The shapefile.prj contains this:
Is there a simple way to find the proj4 code for this?
And I have a netCDF file that I can plot correctly with the following Python code:
map = Basemap(llcrnrlon=lon[0,0], llcrnrlat=lat[0,0], urcrnrlon=lon[560,300], urcrnrlat=lat[560,300], #rsphere=(6378137.00,6356752.3142),\ resolution='l', area_thresh=1E6, projection='lcc', #lat_1=90, lon_0=0) lat_1=71.0, lon_0=-39.0)
I would like to keep the Basemap (python) projection, and convert the shapefile contents to that projection.
Further research and I found this post: Shapefile PRJ to PostGIS SRID lookup table? with a helpful bit of code. Running that and I find the following:
Proj4 is: +proj=stere +lat_0=90 +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs EPSG is: None
I found out I can determine the basemap proj string with:
I'm now using the advice here http://all-geo.org/volcan01010/2012/11/change-coordinates-with-pyproj/ to try to convert from one to the other...