Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  651 1665 (N22E42 L=238)  1578    171        3 ( 576/1616- 737/1701)  33.7 E46-E36 N17-N25
      2 1557 1495 (N35W23 L=303)  1323    251        2 (1441/1435-1660/1546)   2.9 W12-W32 N33-N36
      3  312 1507 (N06E62 L=218)   860    100        1 ( 291/1456- 335/1545) -71.6 E67-E59 N02-N08
      4  441 1432 (N03E46 L=234)   276     42        0 ( 430/1400- 457/1457) -76.0 E48-E45 N01-N05
      5 1813 1204 (N25W48 L=329)  1116    272        1 (1785/1094-1836/1340)  90.0 W45-W51 N17-N33
      6  207 1269 (S09E63 L=217)   642     99        0 ( 188/1226- 235/1316) -32.0 E67-E58 S11-S06
      7  845 1211 (S01E14 L=266)   336     92        0 ( 796/1187- 887/1232)  11.3 E17-E12 S02-N00
      8 1742 1132 (N18W42 L=323)   285     26        0 (1721/1117-1768/1151)  14.0 W40-W44 N16-N19
      9  844  856 (S21E06 L=274)   402    138        0 ( 785/ 815- 899/ 886)  -7.1 E10-E02 S21-S19
     10  508  567 (S47E26 L=254)   805    134        0 ( 450/ 529- 576/ 606)  -4.4 E32-E19 S48-S45
     11  643  531 (S46E13 L=267)   364     48        1 ( 615/ 502- 666/ 561) -63.4 E15-E10 S47-S43
     12 1040  304 (S49W29 L=309)   275     57        0 (1013/ 282-1064/ 331)  -7.1 W24-W32 S50-S48
  

    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