Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  631 1375 (N38E16 L=158)   873    223        1 ( 558/1320- 709/1411)   4.4 E22-E09 N35-N39
      2  162 1208 (N37E61 L=114)   824    270        0 ( 123/1170- 202/1251) -36.9 E64-E56 N33-N39
      3 1798  968 (S19W49 L=223)  1582    285        3 (1757/ 870-1859/1068)  33.7 W43-W57 S24-S14
      4  361 1005 (N21E43 L=132)  1041    198        1 ( 298/ 976- 437/1060) -37.9 E48-E37 N18-N25
      5 1441  787 (S18W18 L=193)   644    177        1 (1374/ 718-1515/ 850)  15.3 W12-W23 S20-S15
      6  198  784 (N10E66 L=109)  1622    265        5 ( 143/ 716- 230/ 845) -76.0 E70-E61 N06-N14
      7  641  562 (S10E35 L=139)  2435    417        2 ( 583/ 385- 734/ 740)  73.1 E39-E29 S23-N00
  

    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