# How to use project to plot the data fall in certain interval of multi cross lines?

**URL:** <https://forum.generic-mapping-tools.org/t/how-to-use-project-to-plot-the-data-fall-in-certain-interval-of-multi-cross-lines/1929>\
**Category:** Q&A\
**Tags:** gmt\
**Created:** [July 26, 2021, 2:28pm UTC](https://forum.generic-mapping-tools.org/t/how-to-use-project-to-plot-the-data-fall-in-certain-interval-of-multi-cross-lines/1929 "2021-07-26T14:28:29Z")\
**Posts on this page:** 8\
**Page:** 1

<div class="post-metadata">

**Author:** ![brioglade](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/brioglade/32/1011_2.png) [@brioglade](https://forum.generic-mapping-tools.org/u/brioglade)\
**Post date:** [July 26, 2021, 2:28pm UTC](https://forum.generic-mapping-tools.org/t/how-to-use-project-to-plot-the-data-fall-in-certain-interval-of-multi-cross-lines/1929/1 "2021-07-26T14:28:29Z")

</div>

This is the script [testtest.txt](https://forum.generic-mapping-tools.org/uploads/short-url/k8x53vB7xmqClEVFpNcip0zQg6S.txt) (1.7 KB)  
This is the data [usgs\_cat2005test.txt](https://forum.generic-mapping-tools.org/uploads/short-url/wW9TdU672iavL8fKh7gehD7Zgus.txt) (3.3 KB)  
This is the result:

 ![profile](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/f/f2e8c798891d436cda29a938fabe05fce9d3c5ed.png)  
The puzzle is that the events projection of the slash cross line does not look correct.  
 ![image](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/f/fe838faa42770f407cb7ac49189f670533bc732b.jpeg)  
i want to get the data projection between the cross lines’ intervals.  
how to do that?Thank you very much.

---

<div class="post-metadata">

**Author:** ![PlanetGus](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/planetgus/32/303_2.png) [@PlanetGus](https://forum.generic-mapping-tools.org/u/PlanetGus)\
**Post date:** [July 26, 2021, 5:49pm UTC](https://forum.generic-mapping-tools.org/t/how-to-use-project-to-plot-the-data-fall-in-certain-interval-of-multi-cross-lines/1929/2 "2021-07-26T17:49:44Z")

</div>

I think what you’re asking for is work in progress (development).

Check this [Grdtrack vertical profile along track](https://flint.soest.hawaii.edu/t/grdtrack-vertical-profile-along-track/1547/)

---

<div class="post-metadata">

**Author:** ![brioglade](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/brioglade/32/1011_2.png) [@brioglade](https://forum.generic-mapping-tools.org/u/brioglade)\
**Post date:** [July 27, 2021, 2:02am UTC](https://forum.generic-mapping-tools.org/t/how-to-use-project-to-plot-the-data-fall-in-certain-interval-of-multi-cross-lines/1929/3 "2021-07-27T02:02:39Z")

</div>

It looks different. What I want is the projection of microseismic event points within a certain range from the cross line.

---

<div class="post-metadata">

**Author:** ![pwessel](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/pwessel/32/30_2.png) [@pwessel](https://forum.generic-mapping-tools.org/u/pwessel)\
**Post date:** [July 28, 2021, 2:01am UTC](https://forum.generic-mapping-tools.org/t/how-to-use-project-to-plot-the-data-fall-in-certain-interval-of-multi-cross-lines/1929/4 "2021-07-28T02:01:56Z")

</div>

I dont understand your drawing. It shows two ~parallel lines while your plot has lines that intersect. Which one is the “slash cross line”?

---

<div class="post-metadata">

**Author:** ![brioglade](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/brioglade/32/1011_2.png) [@brioglade](https://forum.generic-mapping-tools.org/u/brioglade)\
**Post date:** [July 28, 2021, 2:49am UTC](https://forum.generic-mapping-tools.org/t/how-to-use-project-to-plot-the-data-fall-in-certain-interval-of-multi-cross-lines/1929/5 "2021-07-28T02:49:23Z")

</div>

![image](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/b/bc594109fab02609cb61711b3cbae0a4a3dd5752.png)  
I mean the seismic events projection of this cross line may not be correct.It seems that it doesn’t match.  
And what I really want is the projection of microseismic events at a certain distance from both sides or unilateral side of the cross line.Therefore, it can be projected in different continuous areas with the ~parallel cross lines.

---

<div class="post-metadata">

**Author:** ![pwessel](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/pwessel/32/30_2.png) [@pwessel](https://forum.generic-mapping-tools.org/u/pwessel)\
**Post date:** [July 28, 2021, 3:10am UTC](https://forum.generic-mapping-tools.org/t/how-to-use-project-to-plot-the-data-fall-in-certain-interval-of-multi-cross-lines/1929/6 "2021-07-28T03:10:49Z")

</div>

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.

---

<div class="post-metadata">

**Author:** ![brioglade](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/brioglade/32/1011_2.png) [@brioglade](https://forum.generic-mapping-tools.org/u/brioglade)\
**Post date:** [July 28, 2021, 3:31am UTC](https://forum.generic-mapping-tools.org/t/how-to-use-project-to-plot-the-data-fall-in-certain-interval-of-multi-cross-lines/1929/7 "2021-07-28T03:31:16Z")

</div>

![image](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/9/94b5b8cbb5f567ecdfa042584a20ff26f75fcf20.jpeg)

---

<div class="post-metadata">

**Author:** ![pwessel](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/pwessel/32/30_2.png) [@pwessel](https://forum.generic-mapping-tools.org/u/pwessel)\
**Post date:** [July 30, 2021, 3:48am UTC](https://forum.generic-mapping-tools.org/t/how-to-use-project-to-plot-the-data-fall-in-certain-interval-of-multi-cross-lines/1929/8 "2021-07-30T03:48:18Z")

</div>

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_).

 ![one](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/f/f63dda040ea3dacf66856d0014a644cc662bb150.png)