Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  418 1479 (N42E36 L= 99)   421     64        1 ( 389/1444- 447/1506)  45.0 E38-E33 N39-N43
      2  538 1479 (N39E25 L=110)   536    118        1 ( 483/1456- 592/1495) -12.5 E29-E20 N37-N41
      3  962 1314 (N20W02 L=137)   295     58        0 ( 926/1296- 990/1337) -33.7 W00-W04 N18-N20
      4  397 1237 (N26E37 L= 98)   272    108        0 ( 360/1203- 437/1264) -39.8 E40-E34 N23-N28
      5 1700 1074 (S09W44 L=178)   376     68        2 (1678/1035-1725/1106)  36.9 W41-W46 S11-S07
      6  955  808 (S09E08 L=126)  2051    569        2 ( 773/ 618-1218/ 927)   7.9 E22-W07 S17-S03
      7  882  515 (S25E20 L=115)   354     70        0 ( 832/ 502- 920/ 530) -31.0 E22-E17 S26-S23
  

    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