Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1163 1553 (N30W13 L=132)  1365    341        0 (1032/1456-1279/1643) -45.0 W04-W19 N23-N37
      2 1411 1587 (N31W32 L=151)   700    143        0 (1348/1548-1477/1632) -40.6 W28-W36 N28-N35
      3  816 1573 (N34E11 L=108)   238     10        0 ( 796/1554- 835/1591)  45.0 E12-E09 N32-N35
      4  506 1463 (N29E34 L= 85)   264    120        0 ( 458/1432- 562/1488)  20.6 E37-E29 N27-N30
      5  410 1351 (N22E41 L= 78)   933    175        0 ( 376/1301- 472/1413) -53.1 E44-E35 N18-N26
      6  877 1379 (N21E07 L=112)   228     40        1 ( 850/1361- 903/1398)  18.4 E08-E05 N19-N21
      7  164 1348 (N24E71 L= 48)  1678     87        1 ( 134/1300- 192/1394)  63.4 E77-E67 N21-N27
      8  847 1321 (N17E09 L=110)   231     88        1 ( 800/1300- 885/1348)   9.5 E12-E06 N16-N18
      9 1841 1136 (N01W59 L=178)  1237    160        2 (1764/1084-1898/1169) -12.1 W51-W66 S02-N03
     10 1878  877 (S15W64 L=183)  1348    171        1 (1807/ 826-1920/ 937)  21.0 W56-W70 S17-S11
     11  786  849 (S11E16 L=103)   424    120        0 ( 742/ 806- 836/ 896)  39.8 E18-E12 S13-S08
     12 1754  685 (S26W54 L=173)  2360    312        0 (1643/ 579-1885/ 760)  19.8 W44-W70 S33-S21
     13  824  625 (S25E16 L=103)   640    196        1 ( 800/ 560- 886/ 710) -73.3 E17-E12 S29-S19
     14 1560  481 (S40W41 L=159)  1247    167        3 (1494/ 441-1622/ 531)  18.4 W34-W46 S42-S36
     15  501  510 (S30E44 L= 75)   301    114        0 ( 454/ 498- 565/ 523)  -6.3 E47-E38 S31-S29
     16 1376  470 (S39W24 L=142)   662    120        2 (1341/ 414-1404/ 518) -54.5 W20-W27 S43-S35
     17 1388  254 (S58W35 L=154)   484     45        1 (1365/ 238-1418/ 277)  -6.3 W31-W40 S60-S55
  

    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