Timeline for How to reclassify a very large land cover dataset?
Current License: CC BY-SA 3.0
12 events
| when toggle format | what | by | license | comment | |
|---|---|---|---|---|---|
| Jan 7, 2013 at 18:11 | comment | added | DoggoDougal |
@Aaron, whuber made a comment advising against the usage of Reclassify as it might not be the most efficient solution. I have looked into his remark, and found a solution that is very simple and requires no tiling. Thank you for your input, though. I appreciate the help that you provided.
|
|
| Jan 4, 2013 at 17:53 | comment | added | Aaron♦ | @torik Post edited to address your comments. If possible, consider updating your version of ArcGIS to 10.1 SP1. Also, I just noticed that you added "D:\alaska_reclassified.tif" as the outpath for the original script when it should be one file deep and properly formatted, such as: r"D:\temp\alaska_reclassified.tif". | |
| Jan 4, 2013 at 17:52 | history | edited | Aaron♦ | CC BY-SA 3.0 |
Tested the first script with the full dataset and updated the resulting reclassified image.
|
| Jan 4, 2013 at 15:39 | comment | added | DoggoDougal | @Aaron, here is a sample subset, produced using the tiling script that you provided above. The failure occurs when the Reclassify procedure is called. | |
| Jan 4, 2013 at 15:31 | comment | added | DoggoDougal | Here is a snapshot of the Python interpreter's results: pastebin.com/PZXVxae1 | |
| Jan 4, 2013 at 15:24 | comment | added | DoggoDougal | @Aaron, bearing in mind that you provided code to accomplish the tiling, how did you create the subset raster the you used to produce the pictured results? I have completed the SplitRaster tiling (producing 100 subsets of the entire raster dataset), and attempted to loop through them all to reclassify. Reclassification failed, unfortunately, resulting in the same "Unexpected Error" message. | |
| Jan 4, 2013 at 5:04 | comment | added | Aaron♦ | @torik No problem--I'm happy to give my two cents. I think removing the color map is not the way to go. Rather, I would focus on splitting data or 64 bit background processing. | |
| Jan 3, 2013 at 22:13 | comment | added | DoggoDougal | Thank you for the tips, Aaron. I will give it a run as soon as I finish another workaround, which seems to require the removal of the color map (referenced here). This method requires splitting the raster as well, so it makes me wonder if Reclassify original failed due to memory-usage or some other reason. | |
| Jan 3, 2013 at 20:50 | comment | added | Aaron♦ | @DavidF Agreed, there would likely be significant improvement in performance. | |
| Jan 3, 2013 at 20:48 | history | edited | Aaron♦ | CC BY-SA 3.0 |
Added a script to split NLCD dataset prior to reclassification
|
| Jan 3, 2013 at 20:30 | comment | added | DavidF | From a performance standpoint, it would be interesting to try an alternative approach using arcpy.RasterToNumPyArray() and do the reclass in numpy. You would likely want to split the raster up into tiles anyway for memory purposes, but I know that with GDAL, re-classing numpy arrays is very fast. | |
| Jan 3, 2013 at 18:10 | history | answered | Aaron♦ | CC BY-SA 3.0 |