Creating a contiguous raster stack from different sources (vectors and DTM bands)

المشرف العام

Administrator
طاقم الإدارة
This is a conceptual question I've been trying to implement for a while, and not being very good at raster/data processing or python I'd like to hear best practices in dealing with an extent such as the ConUSA.

Problem: I would like to create a tool which would allow me to append csv files with latitude and longitude with values from raster layers and polygons. Basically a spatial join tool but for standard layers...

Ideas: I've been struggling with whether to deal with shapes in postgis or flat files. The problem I have with postgis is I suck at the whole SQLAlchemy engine and such and never know how to properly apply it in a loop when pulling out chunks of shapes by county or state or things like that

My current thought process is this: Take the NLCD2014 in EPSG:5071, and create band with NED 1arc second after doing a gdal_translate of the that gives me 2 bands. Few concerns:

1 - Would it be best to keep the NED in it's 1 degree tiles, project them to 5071 and layer everything into 1 degree x 1 degree tiles. I know it would be easier to process than a merged tif... Are there any good resources or examples of this type of work?

2 - I believe it wouldn't be difficult to import a country wide layer of DFIRM polygons and export as 1 degree by 1 degree tiles... right? Is there a better way?



أكثر...
 
أعلى