Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  810 1852 (N47E46 L=340)  1180    230        1 ( 777/1803- 872/1894) -12.1 E53-E39 N42-N50
      2  601 1450 (N14E36 L=349)  1734    416        3 ( 500/1333- 721/1578) -22.6 E49-E25 N08-N20
      3 1573 1307 (N27W30 L= 55)   628    162        0 (1513/1242-1622/1389) -33.7 W24-W34 N24-N31
      4 1222 1285 (N17W06 L= 32)   429     70        1 (1183/1250-1261/1319) -36.9 W03-W08 N15-N19
      5  605 1149 (S04E28 L=358)   616    126        2 ( 540/1110- 670/1184)  40.6 E31-E24 S07-N00
      6 1241 1129 (N08W10 L= 36)   287     30        0 (1222/1103-1261/1153) -26.6 W08-W12 N07-N10
      7  847  960 (S10E09 L= 17)   732    218        3 ( 730/ 930- 966/1000)  14.9 E16-E01 S11-S07
      8  227  815 (S31E54 L=332)  1110    298        0 ( 171/ 750- 293/ 882)  37.6 E61-E47 S36-S26
      9 1558  839 (S01W37 L= 63)   247    102        0 (1513/ 817-1601/ 866)  -8.1 W33-W40 S01-N00
     10 1300  530 (S25W30 L= 56)   344     63        0 (1266/ 513-1342/ 545)  21.8 W27-W32 S26-S23
     11  881  509 (S37W03 L= 29)   476    143        1 ( 808/ 486- 960/ 534)  24.0 E01-W08 S39-S33
  

    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