Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1089 1382 (N27W07 L=194)   287    100        1 (1052/1354-1138/1403) -26.6 W04-W10 N24-N28
      2  446 1360 (N29E41 L=146)   283     52        1 ( 419/1344- 475/1375)  21.8 E43-E38 N27-N29
      3  561 1340 (N27E31 L=156)   260     69        0 ( 527/1325- 596/1352)  -9.5 E33-E27 N26-N27
      4  359 1092 (N12E46 L=141)   914    187        0 ( 338/1012- 391/1188) -68.2 E48-E42 N07-N17
      5 1661  943 (S06W42 L=229)   952    157        2 (1595/ 902-1739/ 975)   0.0 W36-W49 S07-S03
      6 1851  893 (S11W62 L=249)   488     84        1 (1833/ 860-1876/ 933) -35.5 W59-W66 S13-S08
      7 1559  759 (S16W34 L=221)   241     86        0 (1517/ 730-1585/ 785) -33.7 W30-W36 S18-S13
      8  493  745 (S10E37 L=149)   285     59        1 ( 468/ 724- 519/ 765) -45.0 E39-E35 S11-S08
      9 1658  625 (S26W45 L=231)   725    208        0 (1627/ 539-1697/ 725)  85.6 W42-W46 S32-S19
     10  967  536 (S26E07 L=179)  1996    564        3 ( 833/ 319-1131/ 701)  46.1 E22-W04 S42-S15
     11 1764  604 (S30W58 L=245)   518     63        1 (1726/ 583-1789/ 628)  18.4 W53-W62 S30-S28
     12 1475  579 (S27W29 L=216)   530    123        0 (1436/ 530-1521/ 626)  50.2 W26-W32 S30-S24
  

    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