Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  743 1812 (N42E41 L= 93)  1363    126        3 ( 696/1760- 800/1853) -14.9 E49-E33 N37-N45
      2 1006 1683 (N35E13 L=120)   306     61        1 ( 986/1654-1031/1712) -63.4 E15-E10 N33-N37
      3  487 1644 (N26E51 L= 82)   849    124        1 ( 433/1604- 530/1677) -13.0 E58-E45 N23-N27
      4 1529 1617 (N41W29 L=162)   512    103        1 (1488/1594-1565/1642)  32.0 W24-W33 N38-N44
      5 1602 1419 (N28W33 L=167)   346    114        0 (1565/1368-1633/1461) -45.0 W30-W35 N25-N30
      6 1720 1224 (N18W43 L=177)   688    179        1 (1680/1158-1784/1283) -36.9 W39-W49 N15-N21
      7  131 1000 (S17E70 L= 63)  3851    262        8 (  81/ 903- 199/1083)  70.0 E83-E60 S22-S11
      8  605 1057 (S10E26 L=107)   235     42        0 ( 583/1035- 635/1076)  33.7 E27-E24 S11-S08
      9 1548  998 (N01W32 L=166)   407    140        0 (1492/ 973-1613/1029)   7.1 W28-W36 N00-N03
     10 1859  944 (N04W60 L=194)   513     83        0 (1834/ 903-1881/ 985) -31.0 W56-W63 N02-N06
     11 1486  916 (S05W30 L=163)   227     92        0 (1434/ 905-1530/ 932)   9.5 W26-W32 S05-S03
     12  397  888 (S22E40 L= 93)   354     96        0 ( 357/ 850- 422/ 917) -53.1 E43-E38 S24-S20
     13  876  873 (S17E07 L=127)   246     44        1 ( 844/ 856- 901/ 893)  33.7 E08-E05 S18-S15
     14  168  757 (S32E70 L= 63)   937    103        0 ( 149/ 699- 186/ 806) -90.0 E72-E67 S35-S28
     15 1516  684 (S18W38 L=171)   390     90        1 (1476/ 649-1546/ 722)  80.5 W35-W39 S20-S14
     16 1551  576 (S23W45 L=178)   551     80        1 (1511/ 558-1597/ 600)  18.4 W41-W48 S24-S21
     17  596  427 (S51E27 L=107)   750    123        2 ( 568/ 374- 634/ 477) -58.0 E29-E23 S55-S47
     18 1519  440 (S32W50 L=183)   352     80        0 (1494/ 413-1542/ 464)  63.4 W46-W51 S34-S30
  

    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