У меня есть два шейп-файла, один файл точек и один многоугольник. Они находятся в одной проекции (широта / долгота). Я хотел бы запустить быструю команду из командной строки или сценарий Python с модулем OGR, который создает новый шейп-файл точек только из точек, которые находятся внутри этого многоугольника. Есть ли хороший способ сделать это?