Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  977 1751 (N48W11 L=164)   355     27        0 ( 959/1724- 999/1772)  31.0 W08-W13 N46-N50
      2  857 1568 (N36E03 L=151)   295     77        0 ( 815/1556- 899/1582) -31.0 E05-W00 N34-N37
      3  742 1542 (N36E12 L=142)   364     89        0 ( 709/1503- 776/1579)  38.7 E14-E08 N33-N38
      4 1371  880 (S13W19 L=173)   192     52        0 (1354/ 850-1386/ 907)  56.3 W17-W20 S14-S11
      5 1372  754 (S21W18 L=172)   652    199        0 (1295/ 693-1443/ 832)  28.6 W12-W23 S23-S16
      6 1192  696 (S22W06 L=159)   290    116        0 (1135/ 667-1254/ 723) -35.5 W02-W09 S24-S19
      7  409  645 (S14E47 L=106)   642     71        1 ( 373/ 605- 436/ 688)  45.0 E50-E44 S16-S11
      8  838  590 (S23E18 L=135)   463    178        0 ( 811/ 503- 877/ 662)  48.4 E22-E14 S28-S19
      9  215  554 (S17E77 L= 77)  1902     84        0 ( 191/ 504- 236/ 595)  35.5 E82-E72 S19-S14
     10  542  385 (S32E49 L=104)  1751    251        0 ( 466/ 336- 618/ 444) -24.4 E56-E39 S35-S27
  

    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