Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1368 1630 (N42W34 L=289)   371     74        1 (1327/1615-1405/1649)  -8.1 W30-W37 N40-N43
      2 1238 1607 (N42W21 L=277)   266     48        0 (1214/1585-1255/1630)  31.0 W18-W23 N40-N43
      3 1781 1255 (N13W58 L=314)  1279    180        2 (1725/1175-1838/1312) -38.7 W52-W63 N08-N17
      4  312 1163 (N14E51 L=205)   308     45        0 ( 281/1147- 341/1177)  -8.1 E54-E47 N13-N15
      5 1072  870 (S05W02 L=258)   213    105        0 (1022/ 849-1118/ 885) -21.8 W00-W05 S06-S04
      6  412  847 (S05E42 L=213)   303     58        0 ( 397/ 817- 432/ 882) -76.0 E43-E40 S06-S02
      7 1427  642 (S22W26 L=281)   438     65        1 (1403/ 611-1454/ 671) -68.2 W24-W27 S24-S19
      8  560  560 (S24E35 L=220)   766    203        1 ( 523/ 472- 603/ 625)  65.6 E40-E31 S30-S19
      9 1299  527 (S29W17 L=273)   724    140        2 (1240/ 481-1354/ 575)  32.0 W13-W21 S32-S25
     10  939  499 (S29E08 L=247)   710    271        1 ( 843/ 440-1058/ 553)  23.6 E16-W00 S33-S25
     11 1573  393 (S43W48 L=303)   595     33        0 (1547/ 371-1600/ 415) -21.8 W42-W52 S45-S40
  

    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