Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1688 1097 (N15W43 L=284)  5691   1029       12 (1389/ 783-1882/1413) -27.8 W17-W72 N00-N30
      2  421  923 (S17E38 L=204)   256     61        1 ( 395/ 894- 443/ 950) -56.3 E39-E35 S19-S15
      3  324  893 (S21E46 L=196)   381    113        1 ( 308/ 843- 336/ 943)  90.0 E46-E44 S24-S17
      4 1743  694 (S07W57 L=299)  1231    273        3 (1669/ 663-1830/ 729)  13.2 W50-W68 S10-S04
      5  930  648 (S25W01 L=243)   395    119        1 ( 861/ 624- 985/ 673)   0.0 E03-W05 S26-S24
      6 1461  577 (S20W38 L=279)   546    151        0 (1393/ 556-1518/ 600)   5.7 W32-W42 S20-S18
      7 1385  338 (S37W44 L=286)   504     87        1 (1349/ 313-1414/ 372)  90.0 W41-W46 S39-S34
      8 1302  304 (S41W40 L=282)   328     50        1 (1276/ 289-1324/ 325) -11.3 W38-W43 S42-S39
  

    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