Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  797 1772 (N52E12 L= 69)   913    139        0 ( 734/1727- 875/1815)  28.3 E17-E03 N48-N55
      2  729 1634 (N41E16 L= 64)   953    252        1 ( 632/1606- 828/1667) -16.4 E25-E08 N37-N44
      3  516 1449 (N29E32 L= 49)  1136    260        2 ( 442/1353- 579/1504) -42.5 E38-E26 N22-N33
      4 1096 1332 (N16W07 L= 87)   561    155        3 (1037/1300-1176/1376) -35.5 W03-W11 N14-N19
      5 1001 1299 (N15W01 L= 81)   262     65        1 ( 966/1275-1038/1329) -56.3 E01-W03 N13-N17
      6  237 1282 (N21E54 L= 26)   684    182        3 ( 187/1259- 283/1311) -11.3 E59-E49 N19-N23
      7  124 1013 (N06E69 L= 12)   933    119        1 ( 105/ 970- 152/1050)  29.7 E72-E64 N03-N08
      8  755  796 (S13E18 L= 62)   451     73        2 ( 734/ 751- 778/ 841) -90.0 E19-E16 S15-S09
  

    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