Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  369 1483 (N16E56 L=282)  1647    328        4 ( 276/1380- 451/1627)  72.5 E62-E51 N08-N27
      2  389 1413 (N12E51 L=287)   311     97        1 ( 345/1387- 415/1444)  45.0 E54-E49 N10-N14
      3  159 1174 (S05E70 L=269)  1025    230        1 ( 125/1146- 192/1214) -13.0 E79-E64 S06-S03
      4  776 1148 (N03E17 L=321)   277    110        3 ( 737/1121- 835/1175)   0.0 E19-E12 N01-N04
      5 1320 1104 (N09W17 L=355)   563    162        2 (1273/1052-1387/1175) -24.0 W12-W21 N06-N12
      6 1653 1060 (N12W40 L= 19)   261    101        2 (1640/1007-1665/1114) -66.8 W39-W42 N08-N15
      7  172 1022 (S14E63 L=275)   722     94        3 ( 156/ 983- 192/1071) -68.2 E66-E60 S16-S11
      8 1618  554 (S19W51 L= 29)   490    187        1 (1591/ 509-1656/ 603) -21.8 W46-W56 S21-S15
  

    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