Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1158 1737 (N45W08 L= 98)   292     41        1 (1133/1719-1179/1759)  26.6 W05-W09 N42-N46
      2  470 1721 (N42E56 L= 34)   533     64        2 ( 445/1701- 498/1747) -19.7 E64-E50 N39-N44
      3 1325 1720 (N44W22 L=112)   448     55        0 (1293/1698-1351/1742)  23.2 W18-W25 N42-N46
      4 1395 1477 (N26W23 L=113)   558    176        1 (1334/1416-1461/1538) -41.2 W19-W27 N22-N29
      5  236 1469 (N23E67 L= 23)   563     87        1 ( 216/1427- 260/1505)  78.7 E69-E65 N20-N25
      6 1726 1424 (N25W51 L=141)   370     85        0 (1695/1387-1757/1451) -38.7 W48-W53 N22-N26
      7 1801 1401 (N24W59 L=149)   887    101        2 (1766/1354-1841/1445) -39.8 W55-W63 N21-N27
      8  243 1277 (N10E57 L= 33)   532    105        0 ( 228/1217- 257/1322) -74.1 E59-E55 N06-N13
      9 1844 1074 (N04W58 L=148)  2140    389        4 (1812/ 942-1884/1271) -90.0 W55-W63 S04-N15
     10 1768 1227 (N12W51 L=141)   429     81        2 (1734/1190-1795/1254) -18.4 W47-W53 N10-N14
     11 1713 1154 (N07W46 L=136)   332     64        2 (1689/1126-1734/1187)  63.4 W43-W47 N05-N09
     12 1798 1101 (N05W53 L=143)   512     87        1 (1770/1073-1835/1129)  23.2 W50-W57 N03-N06
     13  783  877 (S14E14 L= 76)   407    160        2 ( 735/ 818- 845/ 913)  40.6 E17-E10 S17-S11
     14  139  838 (S16E71 L= 19)   898    116        1 ( 100/ 800- 185/ 868)   5.4 E84-E63 S17-S14
     15  989  778 (S19E01 L= 89)   405     63        2 ( 968/ 740-1017/ 824)  78.7 E02-W00 S21-S16
     16  946  732 (S22E04 L= 86)   183     61        0 ( 913/ 714- 975/ 749) -11.3 E06-E01 S23-S21
     17  828  669 (S26E12 L= 78)   275     95        0 ( 799/ 636- 866/ 715) -51.3 E13-E09 S28-S23
     18 1136  535 (S34W10 L=100)   352     68        2 (1104/ 516-1174/ 558) -18.4 W07-W13 S35-S32
  

    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