This is the script testtest.txt (1.7 KB)
This is the data usgs_cat2005test.txt (3.3 KB)
This is the result:
I think what you’re asking for is work in progress (development).
Check this Grdtrack vertical profile along track
It looks different. What I want is the projection of microseismic event points within a certain range from the cross line.
I dont understand your drawing. It shows two ~parallel lines while your plot has lines that intersect. Which one is the “slash cross line”?
OK, I will look at the first case to make sure it works or not.
I am still not sure what your other sentences are trying to say - sorry - but we can get to that later.
The output of project -N looks correct to me. I wonder if perhaps you are not using the right coordinates in your projected plots. Given you are using gawk '{if($4<100) print $3,$5,$3}' cross2.dat the coordinates you are using in that plot on the right are the original z-value from USGS column 3 (your “depth”) versus column 5 in cross3.dat, which is the coordinate s from project. The s coordinate is the x-coordinate in the original coordinate system of the projected point after it is moved onto the line q = 0 (the blue circle below). See the figure below. Here the red circle is the original (x,y) coordinates and the blue circle is the projected point whose coordinates are (r,s). The (x,y) point’s projected coordinate (red) is (p,q). So you are plotting the s coordinates of the blue points. Is this what you want? I suspect you meant q. Input (x,y) becomes (p,q) in the projection, while (r,s) is the inverse projection of (p,0).
