Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  858 1822 (N62E23 L= 14)  1192    145        2 ( 793/1797- 928/1849)   0.0 E34-E13 N59-N65
      2 1150 1728 (N53W13 L= 50)   358     43        1 (1123/1707-1175/1745) -11.3 W09-W15 N50-N54
      3  552 1473 (N32E37 L=  0)   491    103        2 ( 518/1430- 589/1521) -36.9 E41-E32 N28-N34
      4 1196 1388 (N27W12 L= 49)   733    144        1 (1132/1338-1241/1423) -29.7 W07-W14 N23-N28
      5  851 1308 (N21E12 L= 26)   285    128        0 ( 808/1279- 903/1339) -18.4 E14-E08 N19-N23
      6 1748  852 (S08W52 L= 89)   529     28        1 (1724/ 826-1773/ 879) -21.8 W49-W54 S09-S06
      7 1699  830 (S10W47 L= 85)   316     76        0 (1689/ 793-1711/ 878) -68.2 W46-W48 S11-S06
      8 1422  756 (S13W26 L= 63)   841    306        0 (1370/ 699-1504/ 835) -77.5 W22-W31 S17-S08
      9  831  687 (S18E12 L= 25)   760    159        2 ( 764/ 654- 906/ 715)  16.7 E17-E07 S19-S15
     10  385  524 (S31E53 L=344)  2664    359        5 ( 304/ 437- 482/ 637) -52.4 E59-E45 S37-S22
     11 1408  571 (S26W27 L= 65)  1153    272        0 (1292/ 537-1531/ 603)   3.0 W18-W37 S28-S23
  

    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