Tag Info

Hot answers tagged

3

A literature search would provide you a wealth of information! Bob McGaughey with USFW-PNW in Seattle, is the developer of FUSION and I am sure would hand over the source code for watershed segmentation. Randy Wynne is at Virgina Tech and developed an IDL virtual machine program implementing a variable window filtering approach Popescu & Wynne (2004). ...


2

Spectral mixture analysis / sub pixel analysis is designed for hyperspectral data, not a 3-band aerial photograph. However, you can try it and see if the output is useful. A tutorial can be found in this pdf and in this ppt/pdf. You will have to skip a significant number of the steps, as you don't have the same amount of information in your dataset. In ...


2

A simple for loop will suffice. You can use readGDAL in the rgdal package but I would recommend raster in the raster package. You have to be a bit tricky and use strsplit in the assign function to strip off the ".tif" file extension. setwd("C:/rasters") rlist=list.files(getwd(), pattern="tif$", full.names=FALSE) for(i in rlist) { ...


1

So thanks to an answer from Mike T to this question, I used listgeo to create TFW files to replicate the geodata that already exists in the TIFs. listgeo -tfw x.tif # FTW Once I did that, MapInfo behaved like a post-Y2K piece of software (although "drag'n'drop" functionality is still beyond it).


1

So GeoTIFFs aren't just GeoTIFF's... Exactly what kind of GeoTIFF have you created? 8bit, 16bit, 24Bit, 32bit? Raster, Grid? ... Sofar MapInfo Pro only supports GeoTIFF Raster up to 24bit, as I recall.


1

I agree with @BradHards comment - it would be best for you to let us know more details. I'll assume you are new to the GIS space, and answer for ease of use. Check out Tilemill from mapbox. It is easy to start with and will allow you to visually set up bounding box, zoom levels, starting zoom, etc. and then you can directly export to .mbtiles


1

A NaN is different than NA. The NaN often results from a divide by zero error whereas NA is the R value for no data. These values behave in specific ways and it would be good for you to read some R background material to understand the behavior. Two useful operators to be aware of are: is.na() and is.nan(). y=c(0,1,2,3,4,NA) x=c(0,1,2,3,4,NA) (d=x/y) ...


1

One strategy is to use a nan_to_num function from Numpy, however it always uses 0. for nan, and there are unfortunately no parameters to replace nans with custom values. $ gdal_calc.py -A nan.tif --outfile=result.tif --calc="nan_to_num(A)"


1

This might be totally tedious but if you convert to Esri ASC (GDAL_Translate -of AAIGRID) which is a text based format and then open in a really good text editor like Notepad++ or VI/VIm, or write a python script, you can replace bad values using find & replace. Then save & close the file and translate to a GeoTiff


1

You can use the flag --NoDataValue=NODATAVALUE to replace NoDataValue. See the third example here


1

If you want to preserve the size of your raster, you have to specify that in the gdalwarp command line. Otherwise gdalwarp tries an own guess, based on the average of x and y source resolution preserving the extent, leading to square raster cells. It is pure coincidence that the pixel and line values are swapped in your case. You can force to preserve the ...


1

Also, you could get the tree canopy volume using this web-based application: LiDARAlphaShape3D http://forest.moscowfsl.wsu.edu:3838/LiDARAlphaShape3D/


1

To visualize, you can check out these web-based applications: LiDARtreesModel3D http://forest.moscowfsl.wsu.edu:3838/csilva/LiDARtreesModel3D/ LiDARstand3D http://forest.moscowfsl.wsu.edu:3838/csilva/LiDARstand3D/


1

You can check out these web-based applicatons: LiDARTreeTop http://forest.moscowfsl.wsu.edu:3838/LiDARTreeTop/ LiDAR3DclusterTree http://forest.moscowfsl.wsu.edu:3838/LiDAR3DclusterTree/



Only top voted, non community-wiki answers of a minimum length are eligible