# GMT5 pscoast issue

**URL:** https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309
**Category:** Q&A
**Tags:** gmt
**Created:** [October 3, 2023, 9:41am UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309 "2023-10-03T09:41:31Z")
**Posts on this page:** 13
**Page:** 1

<div class="post-metadata">

### Author: ![mrav](https://avatars.discourse-cdn.com/v4/letter/m/e79b87/32.png) [@mrav](https://forum.generic-mapping-tools.org/u/mrav)
#### Post date: [October 3, 2023, 9:41am UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/1 "2023-10-03T09:41:31Z")

</div>

I have a problem with pscoast command. I am using the following commands in GMT5 in order to extract a file with the coastlines of Iceland.

**pscoast -R-26/-12/62/68 -Df -W -M \> iceland.txt**

The problem is that I would like to not contain the lakes and the rivers coastlines.  
How can I modify my pscoast command in order to make it??

---

<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: [October 3, 2023, 10:21am UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/2 "2023-10-03T10:21:17Z")

</div>

Does `-A5000` work (plots only area with area greater than 5000m^2)?

---

<div class="post-metadata">

### Author: ![Esteban82](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/esteban82/32/54_2.png) [@Esteban82](https://forum.generic-mapping-tools.org/u/Esteban82)
#### Post date: [October 3, 2023, 3:03pm UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/3 "2023-10-03T15:03:20Z")

</div>

In GMT 6.4 you could use -W1 (I don’t remenber about GMT5).

`pscoast -R-26/-12/62/68 -Df -W1 -M > iceland.txt`

> **[coast — GMT 6.4.0 documentation](https://docs.generic-mapping-tools.org/latest/coast#w)**

---

<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: [October 3, 2023, 4:30pm UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/4 "2023-10-03T16:30:57Z")

</div>

An old-fashion method is to read the man [page](https://docs.generic-mapping-tools.org/6.4/coast.html), e.g. **-A** , and see how you can exclude lakes entirely, limit features (including islands) by area (as @PlanetGus), and more.

---

<div class="post-metadata">

### Author: ![mrav](https://avatars.discourse-cdn.com/v4/letter/m/e79b87/32.png) [@mrav](https://forum.generic-mapping-tools.org/u/mrav)
#### Post date: [October 4, 2023, 5:34am UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/5 "2023-10-04T05:34:44Z")

</div>

Thank you. But I have one more question…Is there a base in which we can see the limit features of all countries ? (except from Iceland) ?

---

<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: [October 4, 2023, 7:15am UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/6 "2023-10-04T07:15:44Z")

</div>

You mean how to not define an area threshold arbitrarily? Not that I’m aware of. But you could probably compute that yourself with binstat (though I’m not sure gmt is the best tool for that)

---

<div class="post-metadata">

### Author: ![mrav](https://avatars.discourse-cdn.com/v4/letter/m/e79b87/32.png) [@mrav](https://forum.generic-mapping-tools.org/u/mrav)
#### Post date: [October 4, 2023, 7:29am UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/7 "2023-10-04T07:29:32Z")

</div>

I mean that for some cases of countries that include plenty of Island (for example Greece) `-A5000` does not working, because I see that Islands coastlines/boundaries are not shown now. Is there a way to not show only the rives and lakes coastlines with pscoast command?

---

<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: [October 4, 2023, 8:44am UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/8 "2023-10-04T08:44:59Z")

</div>

Did you read about -A. It is all there. I am not going to read it for you.

---

<div class="post-metadata">

### Author: ![mrav](https://avatars.discourse-cdn.com/v4/letter/m/e79b87/32.png) [@mrav](https://forum.generic-mapping-tools.org/u/mrav)
#### Post date: [October 9, 2023, 12:40pm UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/9 "2023-10-09T12:40:06Z")

</div>

Thank you but unfortunately in GMT6 it is no use for me. I still see the coastlines of rivers and lakes and I want to ignore them…

---

<div class="post-metadata">

### Author: ![Esteban82](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/esteban82/32/54_2.png) [@Esteban82](https://forum.generic-mapping-tools.org/u/Esteban82)
#### Post date: [October 9, 2023, 12:59pm UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/10 "2023-10-09T12:59:59Z")

</div>

Ok. You should check the documentation of the version that you have.

---

<div class="post-metadata">

### Author: ![mrav](https://avatars.discourse-cdn.com/v4/letter/m/e79b87/32.png) [@mrav](https://forum.generic-mapping-tools.org/u/mrav)
#### Post date: [October 9, 2023, 1:11pm UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/11 "2023-10-09T13:11:20Z")

</div>

I mean I have downloaded the GMT6 but still no use… Thank you

---

<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: [October 9, 2023, 1:23pm UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/12 "2023-10-09T13:23:40Z")

</div>

Clearly you’re not trying to read the documentation.  
`-A` is exactly what you are asking (`-A0/0/1`)

---

<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: [October 9, 2023, 1:24pm UTC](https://forum.generic-mapping-tools.org/t/gmt5-pscoast-issue/4309/13 "2023-10-09T13:24:04Z")

</div>

Are you telling us that despite using **-A0/0/1** you are seing lakes. This was implemented in GMT 3 so I doubt your statement.
