Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1324 1852 (N54W40 L= 67)   490    167        1 (1297/1820-1354/1879)  26.6 W34-W46 N50-N56
      2 1041 1634 (N37W06 L= 33)   262     49        0 (1023/1611-1061/1656) -90.0 W04-W07 N35-N38
      3 1769 1425 (N18W57 L= 84)  2500    473       11 (1706/1359-1839/1490) -35.0 W50-W64 N14-N22
      4  996 1399 (N21W01 L= 28)   207     42        0 ( 969/1381-1024/1422) -56.3 W00-W02 N19-N22
      5 1236 1219 (N08W14 L= 41)   202     44        0 (1202/1207-1266/1232) -18.4 W11-W15 N07-N09
      6  432 1045 (N03E38 L=350)  1126    252        1 ( 374/ 920- 470/1148) -74.1 E41-E35 S04-N10
      7  251 1073 (N07E53 L=334)   352     95        2 ( 227/1035- 290/1111)  45.0 E55-E49 N04-N08
      8  189  999 (N03E60 L=327)  2721    539        8 ( 126/ 921- 284/1081) -23.6 E68-E50 S02-N07
      9  279  773 (S11E54 L=333)   751    191        2 ( 232/ 725- 324/ 831)  24.4 E59-E48 S13-S08
     10  307  590 (S23E56 L=331)   752    137        2 ( 265/ 558- 350/ 628) -33.7 E60-E52 S25-S20
     11  982  402 (S42E09 L= 18)   795    192        1 ( 900/ 360-1084/ 460)  18.4 E15-W00 S44-S37
     12  794  344 (S44E26 L=  2)   299     69        1 ( 759/ 331- 830/ 363)  -9.5 E28-E22 S45-S42
  

    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