Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  781 1648 (N27E29 L=195)   713    209        0 ( 725/1561- 817/1741) -38.2 E37-E23 N21-N32
      2  472 1492 (N12E45 L=179)   345     74        1 ( 437/1470- 504/1508)  11.3 E48-E42 N10-N13
      3  262 1292 (S02E57 L=167)  1625    415        2 ( 237/1110- 301/1485) -61.4 E64-E49 S12-N10
      4 1370 1304 (N16W15 L=239)   800    326        3 (1190/1276-1498/1326)  24.0 W04-W23 N10-N19
      5  382 1038 (S16E42 L=182)   783    287        1 ( 284/ 951- 461/1146)  45.0 E49-E36 S21-S08
      6  225 1091 (S14E56 L=167)   583     93        4 ( 194/1054- 252/1135) -45.0 E60-E53 S15-S11
      7  904 1023 (S09E07 L=217)   229     51        1 ( 879/1000- 927/1049)  33.7 E08-E05 S10-S07
      8 1207  856 (S13W14 L=238)   301    126        0 (1158/ 838-1259/ 886)   9.5 W10-W16 S14-S10
      9  324  726 (S35E47 L=177)   875    130        1 ( 278/ 681- 386/ 764)  24.8 E53-E40 S38-S32
     10 1634  479 (S24W58 L=281)   341     68        0 (1608/ 444-1656/ 507)  90.0 W56-W58 S26-S22
  

    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