Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1561 1708 (N56W52 L=159)   843    106        3 (1531/1671-1591/1732)  20.2 W45-W65 N52-N59
      2 1010 1393 (N24E05 L=102)   862    409        0 ( 919/1275-1132/1502) -39.1 E13-W04 N17-N30
      3  348 1453 (N20E58 L= 49)   444    101        2 ( 324/1425- 380/1497)  78.7 E59-E54 N17-N22
      4 1514 1426 (N32W32 L=138)   495    125        0 (1460/1377-1548/1467) -39.8 W27-W34 N28-N34
      5  920 1313 (N18E10 L= 97)   378    150        2 ( 906/1230- 939/1371) -77.5 E11-E08 N12-N21
      6 1134 1312 (N20W04 L=111)   648    172        1 (1069/1248-1204/1341) -20.6 W00-W08 N16-N21
      7  841  760 (S17E09 L= 98)   520    116        1 ( 792/ 727- 902/ 790)  29.7 E11-E04 S19-S14
      8 1390  608 (S21W30 L=137)   387    130        2 (1365/ 572-1423/ 649) -68.2 W27-W32 S23-S18
      9  256  594 (S36E65 L= 41)   701     89        0 ( 223/ 563- 296/ 618)  12.8 E79-E57 S39-S33
     10 1440  386 (S36W44 L=150)   715    132        1 (1393/ 345-1492/ 427)  66.8 W40-W46 S39-S32
  

    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