Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  807 1735 (N34E33 L= 94)   982    142        1 ( 777/1663- 834/1801) -39.8 E39-E27 N28-N38
      2  500 1744 (N30E61 L= 65)  1064    114        2 ( 465/1692- 539/1787) -17.7 E73-E51 N26-N33
      3  676 1606 (N22E35 L= 91)   447    130        1 ( 620/1558- 734/1644)  56.3 E37-E32 N18-N25
      4 1227 1617 (N34W01 L=127)   317     80        0 (1200/1595-1260/1643) -26.6 E01-W03 N32-N35
      5  391 1572 (N16E57 L= 69)   491    112        1 ( 366/1547- 429/1612)  -7.1 E60-E52 N14-N19
      6  993 1472 (N19E11 L=115)  1613    482        5 ( 839/1377-1144/1605) -14.6 E23-W00 N14-N25
      7  708 1377 (N08E26 L=101)  1328    383        1 ( 676/1221- 734/1582) -67.8 E31-E21 S01-N21
      8 1514 1438 (N29W24 L=150)   800    253        0 (1468/1378-1604/1494)  35.0 W20-W30 N23-N32
      9 1635 1377 (N28W34 L=160)   410    131        0 (1608/1315-1660/1432) -56.3 W31-W36 N24-N30
     10 1015  914 (S13W02 L=128)   684    202        2 ( 937/ 886-1120/ 948)  29.1 E02-W07 S16-S10
     11  871  851 (S20E06 L=120)   236     30        2 ( 848/ 824- 895/ 877) -71.6 E07-E04 S21-S18
     12  592  832 (S26E24 L=102)   364    126        0 ( 531/ 806- 636/ 854)   0.0 E28-E20 S26-S24
     13  555  694 (S35E26 L=100)   361     80        0 ( 525/ 674- 597/ 722)  -9.5 E28-E22 S36-S33
     14 1412  499 (S29W39 L=165)  1104    246        3 (1326/ 408-1477/ 578) -16.7 W29-W49 S32-S25
     15 1223  529 (S31W25 L=151)   340     91        0 (1168/ 515-1265/ 543)  15.9 W20-W27 S32-S29
     16  719  489 (S46E11 L=116)   717    104        2 ( 669/ 467- 788/ 521)  15.3 E15-E04 S47-S43
     17 1410  373 (S36W47 L=174)   541     72        1 (1381/ 341-1451/ 405)  78.7 W44-W49 S38-S33
     18  738  237 (S68E04 L=122)   971    137        3 ( 673/ 210- 799/ 265)  -2.9 E15-W06 S69-S65
  

    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