Tôi đang cố gắng tính toán hồ sơ Sérsic của các thiên hà khác nhau từ SDSS dựa trên các hình ảnh được cung cấp bởi trang web của sở thú thiên hà. Tôi đang làm điều này như là một phần của cuộc thi kaggle về việc sử dụng máy học để dự đoán hình thái thiên hà. Tôi không có cơ hội xếp hạng cao trong cuộc thi này vì vậy tôi không ngần ngại yêu cầu giúp đỡ.
Tôi đã sử dụng hàm R contourLines để xác định các isophot của thiên hà và sau đó khớp các hình elip cho mỗi isophote. Điều này dường như hoạt động tốt, các isophotes hầu như luôn luôn phù hợp bởi các hình elip và các hình elip gần như đồng tâm. Sau đó, hãy để Icường độ điểm ảnh của một isophote và Rlà độ dài của trục bán chính của hình elip tương ứng, tôi cần phải phù hợp với một phương trình có dạng
log I(R) = log I_0 - k * R^(1/n)
Cách tiếp cận đơn giản dường như là lấy nhật ký của cả hai bên và sử dụng hồi quy OLS, vì vậy tôi phù hợp với mô hình tuyến tính trong R của mẫu
log(log(I)) ~ log(R)
Các biểu đồ kết quả cho thấy sự phù hợp tốt nhưng các chỉ số Sérsic n kết quả hầu như luôn nhỏ hơn một và không bao giờ lớn bằng hai. Điều này có vẻ không đúng vì các chỉ số từ 4 trở lên dường như phổ biến trong cách đọc của tôi. Tôi không nhận được bất cứ nơi nào gần 4 cho hình ảnh của M87.
Có thể lấy nhật ký nhật ký làm phẳng mọi thứ quá nhiều và chỉ mục không đủ đáp ứng. Tôi đã thử sử dụng nls để làm việc chỉ với nhật ký nhưng nó không di chuyển các chỉ số nhiều.
Có phần mềm hay thuật toán tiêu chuẩn nào để tính toán chỉ số Sérsic từ một hình ảnh không? Có hình ảnh tham khảo nào tôi có thể làm việc từ đó sẽ cho phép tôi kiểm tra xem thuật toán của tôi có hợp lý không? Bất kỳ khuyến nghị về cách tiến hành sẽ được hoan nghênh.
CẬP NHẬT: Tôi đã tìm thấy các chương trình GALFIT và GIM2D trông giống như chúng có thể hữu ích. Bất kỳ phần mềm khác thường được sử dụng cho việc này?