Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  305 1535 (N06E68 L=246)  1475    135        2 ( 259/1479- 348/1606)  71.6 E71-E65 N02-N11
      2  916 1419 (N14E17 L=296)   254     67        0 ( 889/1383- 941/1450)  78.7 E18-E16 N10-N15
      3 1628 1109 (N17W34 L=347)   391    103        1 (1602/1073-1665/1137)  59.0 W31-W36 N14-N19
      4 1661  903 (N07W42 L=355)  1895    344        6 (1567/ 807-1722/ 984) -21.0 W33-W47 N02-N11
      5 1830  813 (N08W60 L= 14)   918    194        2 (1801/ 738-1852/ 881) -28.6 W54-W66 N04-N10
      6  531  720 (S36E23 L=290)   460    139        1 ( 478/ 688- 592/ 743)  36.9 E27-E19 S39-S33
      7 1666  443 (S17W65 L= 18)  1135    229        1 (1642/ 393-1688/ 503) -24.4 W60-W71 S19-S13
  

    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