GIS Rasterize Vector Data "requires numeric data"
Hello,
I am trying to take a HYSPLIT particle concentration map and import it into Igor to show contour lines of equal concentration over a region. HYSPLIT has an add on to export data as a .shp file (not its native file type), so I have been using that file type as a go between for the two softwares. The .shp file HYSPLIT generates is easily loaded in with the GISLoadVectorData command as detailed in the igor help topic, but all of the waves that are created (except W_geometries and W_geometriesIndices) are imported as text waves. When I try to use GISRasterizeVectorData, it alerts me that "Rasterizing vector data requires numeric data." Even when I convert the text wave to a numeric wave using the function below, I get the same error. I believe rasterizing this data is the best way to display the data, does anyone have suggestions to get this to work? Attached is my working file (the forum won't allow me to include the .shp file).
Thank you in advance,
Walt
Function ConvertTxtWvToNumWv(W)
Wave /T W
Variable np = numpnts(W)
Make /O /N=(np) NumConvert
Variable x
for(x=0;x<np;x+=1)
NumConvert[x] = str2num(W[x])
endfor
End
Forum
Support
Gallery
Igor Pro 9
Learn More
Igor XOP Toolkit
Learn More
Igor NIDAQ Tools MX
Learn More
What is your goal? It looks like what you loaded contains traces for a contour plot. You have already superposed those contours on top of a map of the US.
Unfortunately, the fellow that wrote the GIS package is no longer with WaveMetrics, so we may have to learn along with you :)
May 30, 2019 at 09:05 am - Permalink
Is the .shp file too large to post? You could share it with WaveMetrics support via DropBox or similar. If we can't figure out the answer to your problems here, you may need to move this to support@wavemetrics.com.
May 30, 2019 at 09:06 am - Permalink
My motivation for the raster is so I can assign colors and labels to each contour individually, creating a display similar to the last image in the DisplayHelpTopic "IgorGIS Package" > IgorGIS Tutorial: Rasterizing Vector Data (the example about geography of the US). Eventually, I will need to create several layers of concentration data for different altitudes at the same XY coordinates.
The .shp file is small (12 KB), but the forum doesn't support the file type. I can email it or transfer it by other means if necessary.
May 30, 2019 at 09:43 am - Permalink
I take it that the vector data you want is the same as that which made the horseshoe-shaped contours. That pretty has to be numeric data, so I'm not sure what the problem is. Since GISRasterizeVectorData requires a fieldname input, is there a chance that you used to wrong name there? I suppose there may also be the possibility that the error message is simply not being helpful. I presume you used the same GISFileHandle as you did to load the data originally.
Another possibility- if the raw concentration data is available, you could make the contour plot with Igor instead of by loading the contours from the .shp file. Is that at all possible?
May 30, 2019 at 01:29 pm - Permalink
Sorry for the delayed response. The horseshoe shaped contours are the vector data I want rasterized. They are easy to display by plotting W_geometries[1] vs W_Geometries[0], but then I have no way to display which contour band corresponds to which concentration. I will look further into raw HYSPLIT concentration data, but the output formats are very limited.
I am trying to use the following command. The issue is that the CONC wave is listed as a string in W_AttributeNames, even though all the values are integers. I am using the same GISFileHandle (ConcHandle) as I did to load the data originally.
GISRasterizeVectorData /D=root:HYSPLIT:fireData ConcHandle, "CONC"
June 3, 2019 at 06:37 am - Permalink
Let's do this as a tech support incident- send an e-mail to support@wavemetrics.com and attach a copy of the .shp file. In that message, it would be handy if you could repeat some of what you posted here, especially including the GIS commands you have used to load the data and that you attempted to use for rasterizing.
June 3, 2019 at 09:48 am - Permalink