Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1539 1683 (N35W49 L=180)   403     77        0 (1515/1661-1560/1705)  23.2 W45-W52 N33-N36
      2  521 1513 (N35E31 L=100)  2489    672       11 ( 311/1431- 724/1602) -23.0 E55-E15 N27-N45
      3  280 1193 (N17E49 L= 82)   793    285        3 ( 257/1132- 331/1249)  45.0 E51-E44 N13-N20
      4  903 1174 (N09E06 L=126)   220     86        0 ( 866/1136- 943/1205) -51.3 E07-E03 N06-N11
      5 1830 1022 (S09W56 L=188)   644    174        4 (1795/ 971-1867/1053)  33.7 W52-W60 S11-S06
      6 1830  902 (S16W56 L=187)  1223    311        4 (1747/ 822-1903/ 984)  25.2 W47-W64 S20-S11
      7  952  788 (S14E07 L=124)   222    110        0 ( 910/ 743-1015/ 821)  23.2 E09-E02 S16-S12
      8  920  471 (S34E14 L=117)   530    220        1 ( 880/ 403- 948/ 581)  59.7 E18-E10 S38-S26
      9  606  378 (S37E41 L= 90)  1099    413        1 ( 581/ 245- 651/ 524)  60.3 E47-E34 S47-S26
     10 1291  452 (S40W13 L=144)   647    149        3 (1210/ 417-1363/ 491)   9.5 W06-W18 S42-S37
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y