Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  543 1636 (N59E20 L= 97)   382     64        0 ( 519/1614- 582/1653)   0.0 E22-E14 N56-N60
      2 1167 1623 (N35W30 L=147)  1152    209        2 (1075/1598-1275/1652) -30.3 W23-W36 N31-N38
      3  287 1348 (N44E50 L= 67)   545     87        0 ( 255/1314- 321/1371) -28.6 E55-E44 N40-N46
      4 1749 1179 (S07W50 L=167)   322     52        1 (1725/1149-1765/1204) -76.0 W47-W51 S08-S04
      5 1402 1042 (S04W21 L=138)   582    148        0 (1373/ 975-1439/1115)  82.9 W20-W22 S08-N00
      6 1718  979 (S17W42 L=159)  1192    305        1 (1625/ 884-1810/1080)  14.0 W32-W52 S19-S14
      7  248  800 (N11E59 L= 58)   844    132        0 ( 221/ 737- 271/ 854)  74.1 E61-E55 N07-N14
      8 1463  554 (S35W14 L=131)  1203    232        0 (1406/ 470-1495/ 656)  67.4 W08-W17 S40-S27
      9  654  434 (S19E41 L= 76)   410    102        0 ( 637/ 391- 672/ 493)  54.5 E44-E38 S21-S14
     10  592  359 (S22E50 L= 67)   708    105        0 ( 550/ 318- 635/ 391)   5.2 E57-E46 S24-S20
  

    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