Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  660 1714 (N24E47 L=240)   730     97        1 ( 615/1683- 708/1737)   0.0 E50-E41 N22-N25
      2 1340 1615 (N38W01 L=288)  1424    431        1 (1168/1515-1514/1703)  -2.0 E13-W16 N35-N39
      3  764 1451 (N10E27 L=260)   784    244        1 ( 686/1368- 828/1550) -26.6 E34-E20 N07-N14
      4  361 1488 (N04E57 L=229)   888    124        2 ( 324/1433- 386/1547) -54.5 E63-E52 N00-N07
      5 1598 1253 (N23W28 L=314)   335    114        0 (1550/1226-1641/1280)   0.0 W24-W31 N21-N23
      6  127 1121 (S21E72 L=214)   680    141        1 ( 103/1052- 163/1174)  45.0 E76-E66 S25-S17
      7 1832 1136 (N25W51 L=337)   553    136        0 (1784/1087-1868/1160)  -7.1 W45-W55 N22-N26
      8 1868  889 (N12W61 L=347)  3007    381        7 (1774/ 796-1944/1034) -15.1 W47-W79 N08-N17
      9  853  960 (S14E08 L=279)   275     44        1 ( 829/ 934- 883/ 986)  90.0 E08-E05 S16-S12
     10 1214  810 (S13W16 L=303)   440    132        0 (1173/ 757-1256/ 873) -20.6 W12-W20 S14-S10
     11 1118  727 (S20W14 L=300)   345    127        0 (1052/ 709-1182/ 743)  14.0 W09-W17 S21-S18
     12 1834  712 (N02W65 L=351)   757     78        1 (1787/ 691-1874/ 732)  20.0 W59-W70 N00-N04
     13  266  659 (S48E54 L=232)   492     65        1 ( 235/ 645- 301/ 678)   9.5 E60-E48 S49-S46
     14  986  652 (S28W09 L=295)   482    204        0 ( 912/ 613-1076/ 672)  26.6 W05-W15 S31-S25
     15 1592  528 (S16W52 L=339)   882    164        3 (1520/ 506-1680/ 560)  26.6 W46-W59 S18-S12
     16  708  467 (S48E03 L=283)   385     19        1 ( 681/ 443- 732/ 492) -21.8 E05-W00 S49-S46
     17  837  440 (S46W08 L=294)   309     81        0 ( 796/ 421- 887/ 459)   7.1 W04-W12 S47-S44
     18  990  379 (S44W22 L=308)   482     99        1 ( 956/ 350-1039/ 418) -12.5 W17-W26 S45-S42
     19 1460  323 (S30W59 L=345)  1157    202        2 (1365/ 273-1534/ 370)   8.5 W48-W70 S33-S28
  

    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