How can I convert XYZ format data to create a line of sight map?

I have a problem determining line of sight using LOS2 or VIEWSHED. I have DEM data with vectors each (of equal lengths) for latitude, longitude and height. Now I want to create a line of sight map and believe that VIEWSHED or LOS2 will help. However, they require map data in terms of Z and R. How can I convert data such that I may be able to use these functions?

Answer 1

[Chaitanya Chitale]

The data are scattered, that is, points here and there in latitude and longitude. The data have to be gridded before can use functions like viewshed, los2, contourfm, etc. We do not yet have a function to do this. However, there is an example that shows the kind of processing that is needed. See the following code:

load seamount
lon = x;
lat = y;
delta = 0.01;
loni = (210.8 : delta : 211.8);
lati = (-48.5 : delta : -47.9)';
Zi = griddata(lon, lat, z, loni, lati);
R = makerefmat(loni(1), lati(1), delta, delta);
figure('Color','white')
worldmap(Zi, R)
geoshow(Zi, R, 'DisplayType', 'surface')