Спасибо тебе за ответ! Я уже нашел решение раньше, поэтому выложу его здесь для тех, кому интересно:
сначала я создаю достаточное количество сэмплов пуассоновского диска в единицах квадратов (достаточно означает больше, чем n)
затем я сортирую эти выборки по меньшей координате (например, точка (10,9), меньшая координата равна 9 - другая точка (8,50), меньшая координата равна 8 - порядок точек будет (8,50), (10,9))
тогда я беру первые n образцов в отсортированном списке. из-за режима сортировки эти образцы снова будут лежать в квадратной области. Затем я масштабирую координаты так, чтобы они снова лежали в единичном квадрате. Теперь у меня ровно n образцов пуассоновских дисков на единицу площади.
затем я использую отображение плоскости на сферу, описанное в http://www.cs.rutgers.edu/~decarlo/readings/mcrt-sg03c.pdf стр. 23, чтобы получить равномерно распределенные выборки на сегменте сферы произвольного угла области
хорошо работает для меня