Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  700 1813 (N39E47 L=329)   949     94        2 ( 660/1789- 749/1840)  -8.7 E54-E41 N36-N41
      2  950 1746 (N38E22 L=354)  1417    196        6 ( 864/1678- 999/1789) -21.8 E31-E16 N33-N41
      3  354 1458 (N09E55 L=322)   742    155        1 ( 294/1413- 409/1501)  32.0 E59-E50 N06-N12
      4  665 1439 (N12E30 L=346)   344     69        1 ( 637/1408- 698/1470) -21.8 E32-E27 N10-N13
      5  893 1411 (N14E15 L=  1)   311     98        1 ( 847/1382- 931/1429)  -9.5 E18-E12 N12-N14
      6 1748 1223 (N20W45 L= 61)   324     61        0 (1715/1205-1780/1236)   8.1 W41-W48 N19-N21
      7  404 1016 (S17E40 L=336)   296     42        1 ( 392/ 986- 422/1050) -76.0 E41-E38 S18-S14
      8  174  783 (S33E68 L=309)  1117    121        1 ( 145/ 725- 198/ 843)  25.0 E77-E62 S36-S29
      9 1120  692 (S24W13 L= 29)   705    206        1 (1097/ 591-1152/ 794) -72.9 W10-W15 S30-S17
     10 1369  528 (S28W35 L= 51)  1197    351        1 (1245/ 458-1507/ 643)  55.5 W28-W39 S34-S18
     11 1569  520 (S24W50 L= 67)   999    306        0 (1455/ 435-1683/ 603)  56.3 W45-W55 S31-S16
  

    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