Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  410 1090 (N24E39 L=318)   598    172        2 ( 348/1056- 486/1120) -31.0 E43-E32 N20-N26
      2  328  974 (N18E47 L=309)   735    145        0 ( 282/ 935- 387/1018) -36.9 E50-E42 N15-N21
      3 1413  927 (S08W20 L= 17)   431    120        0 (1362/ 887-1469/ 952)   0.0 W16-W23 S09-S06
      4  692  811 (N02E24 L=332)   818    189        1 ( 667/ 721- 721/ 899)  71.6 E27-E21 S02-N07
      5 1238  784 (S12W07 L=  3)   724    171        2 (1199/ 692-1267/ 872)  59.0 W02-W09 S16-S06
      6 1580  776 (S21W29 L= 25)   502    172        0 (1554/ 696-1594/ 867)  60.9 W25-W30 S25-S16
      7 1350  770 (S15W13 L=  9)   231     47        0 (1333/ 740-1371/ 798)  45.0 W11-W14 S16-S13
      8 1631  631 (S33W31 L= 27)   308     53        0 (1598/ 602-1652/ 655)  11.3 W28-W33 S33-S31
      9 1217  528 (S27E01 L=356)   870    244        0 (1099/ 436-1330/ 621)  18.4 E10-W08 S30-S23
     10  877  448 (S24E24 L=332)   342     91        0 ( 841/ 414- 910/ 487)  23.2 E27-E20 S25-S21
     11  851  225 (S40E39 L=318)   722    178        1 ( 832/ 171- 872/ 267)  41.6 E43-E33 S45-S36
  

    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