I currently used the knickpointfinder MATLAB codes from TopoToolbox. I want to plot the locations of knickpoints (which are shown over longitudinal profile) over DEM. Please guide me with the MATLAB codes for plotting knickpoint position over DEM.
I was dealing with the same issue,
It is easier when only using the main trunk knickpoints...
The code should be something like this, where points are the knickpoints coordinates x (distance along profile) and y(elevation), CAT is a structure where I stored several stream objects, level is each of the profiles. Some variables should be modified is just a draft:
[row,col]=find(CAT.STREAMOBJ.distance==points(1));
Xknk(level)=CAT.STREAMOBJ.x(row,1);
Yknk(level)=CAT.STREAMOBJ.y(row,1);
Zknk(level)=points(2);
levelo(level)=level;
%here export the kniuckpoints to shapefile to display in arches
[Tracks(1:length(Xknk)).Geometry] = deal('Point');
for i = 1:length(Xknk)
[Tracks(i).Catchment]= Catchment;
[Tracks(i).level]= levelo(i);
[Tracks(i).X]= Xknk(i);
[Tracks(i).Y]= Yknk(i);
[Tracks(i).Z]= Zknk(i);
end
fout=('Knickpoints.shp');
shapewrite(Tracks, fout);
Hello.
You could use the ArcGis profiler toolbar extension, which is recognized by Matlab.
http://geomorphtools.geology.isu.edu/Tools/StPro/Tutorials/StPro_UserGuides_Final.pdf
Kind Regards.
Hello.
You can plot the knickpoints on map view using the "kp" variable created with the knickpointfinder script of topotoolbox.
this is the code i use:
%%Knickpoint map
figure(4)
imageschs(DEM,[],'colormap',[.9 .9 .9],'colorbar',false);%hillshade
hold on
plot(S, 'b', 'LineWidth',1) %plots the stream netwkork
scatter(kp.x,kp.y,50,'MarkerFaceColor','r'); %plot the knickpoints X and Y coordinates
hold off
title({['Knickpoint and streams map ']})%%Título en dos lineas
i hope it will be usefull for you.
- Badges
- Science topic