Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1569 1588 (N20W54 L= 50)  1343    190        2 (1493/1549-1646/1615) -35.0 W48-W59 N15-N22
      2 1717 1249 (S02W50 L= 46)   838    128        1 (1689/1194-1747/1325)  63.4 W47-W52 S05-N01
      3  708 1020 (N14E18 L=338)   233     93        0 ( 667/ 999- 754/1035)   0.0 E20-E14 N13-N14
      4  219  880 (N16E59 L=298)   350     53        1 ( 202/ 849- 239/ 899)  71.6 E60-E56 N14-N17
      5  265  761 (N09E58 L=298)  1693    398        5 ( 149/ 708- 361/ 811) -17.7 E72-E50 N04-N12
      6  934  523 (S20E20 L=337)   395    131        0 ( 903/ 466- 959/ 582)  35.5 E23-E16 S23-S17
      7  541  313 (S24E57 L=300)  2194    312        3 ( 500/ 224- 587/ 406)  54.5 E61-E50 S31-S17
      8 1050  248 (S43E27 L=330)   360     72        0 (1030/ 219-1068/ 275)  33.7 E29-E23 S45-S41
  

    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