Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1310 1391 (N25W07 L= 36)   791    232        2 (1251/1279-1375/1486) -37.6 W01-W14 N19-N29
      2  856 1300 (N07E17 L= 12)   233     64        0 ( 823/1285- 893/1313)  36.9 E18-E14 N05-N08
      3  498  894 (S24E30 L=359)   861    225        1 ( 411/ 859- 583/ 948)  33.7 E36-E24 S27-S19
      4 1133  832 (S12W11 L= 40)   268     99        2 (1091/ 814-1189/ 854)  45.0 W09-W13 S13-S09
      5 1182  786 (S13W15 L= 44)   217     86        0 (1141/ 775-1231/ 799)  36.9 W12-W17 S14-S11
      6  537  643 (S39E23 L=  6)  1496    244        2 ( 455/ 549- 627/ 738) -23.6 E30-E14 S43-S34
      7  898  711 (S25W01 L= 30)   342     96        0 ( 845/ 689- 954/ 738)   9.5 E02-W04 S25-S24
      8 1599  693 (S07W45 L= 74)   499     97        1 (1543/ 675-1634/ 719)  33.7 W41-W47 S09-S04
      9 1465  608 (S15W39 L= 68)   473    136        0 (1433/ 545-1505/ 673) -69.4 W37-W41 S18-S10
     10  787  562 (S38E02 L= 27)   424     85        0 ( 740/ 544- 837/ 580)   8.1 E05-W02 S39-S36
     11 1319  540 (S23W32 L= 61)   529    244        0 (1214/ 523-1431/ 553)  18.4 W25-W40 S26-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