HI all,
Of course I am going to offer up a GHS solution..... but how you do it will depend on whether you want high-contrast stars (with a Gaussian shape), or low contrast shape (round with little/no halos).
Using GHS with a high b (maybe a couple of stretches with HP to control brightness), will give the former, as suggested above.

For the latter, check the invert box in GHS, apply an intermediate +ve b, move SP (right up to 1 initially), and apply stretch to form circular stars. To lower b towards 0 and adjust D to achieve your desired star shape. This will leave the star brightness nearly saturated (so they won't really hold colour saturation). To reduce star brightness, reduce SP from 1 to get the desired star brightness. With the controls you can somewhat independently control the star brightness, size, and colour retention.


I often use the latter technique to create a mask from my extracted NB star_mask luminance to past RGB stars into my starless NB image. Alternatively you can blend a starry version onto a starless version of luminance image with the pixelmath: Starryimage = Starlessimage + (1-Starlessimage)*star_mask and then LRGB combine with your starry colour version.
CS & I hope this makes sense...
Dave