Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1123 1831 (N47E23 L=110)   366     41        0 (1099/1811-1147/1848) -11.3 E26-E20 N45-N49
      2  947 1799 (N38E32 L=101)  1361    162        0 ( 884/1762-1030/1826)  68.2 E36-E27 N33-N43
      3  540 1701 (N21E55 L= 79)   754    101        2 ( 508/1657- 580/1736) -68.2 E58-E50 N18-N24
      4 1017 1593 (N26E17 L=117)   318     61        0 (1004/1557-1031/1630) -51.3 E18-E14 N23-N28
      5 1370 1546 (N33W06 L=140)   246     29        1 (1354/1524-1386/1568) -33.7 W04-W07 N31-N34
      6 1564 1378 (N29W24 L=157)   608    153        0 (1539/1312-1597/1453) -58.0 W20-W27 N25-N33
      7 1276 1390 (N20W05 L=138)   198     78        0 (1256/1349-1296/1424) -36.9 W02-W06 N18-N22
      8  817 1391 (N08E21 L=112)   691    179        0 ( 731/1373- 909/1411)  21.8 E26-E16 N06-N10
      9  264 1148 (S17E54 L= 79)  1341    231        2 ( 218/1055- 309/1236)  71.6 E57-E50 S23-S11
     10  758 1152 (S07E18 L=115)   415    111        1 ( 720/1105- 796/1200) -33.7 E20-E14 S08-S04
     11  529 1162 (S11E33 L=101)   452     92        0 ( 500/1125- 567/1195)  59.0 E34-E30 S13-S08
     12 1688 1011 (N11W40 L=174)   344     49        1 (1658/ 985-1713/1040)  53.1 W38-W41 N09-N13
     13  801  984 (S15E12 L=122)   798    276        1 ( 724/ 934- 925/1037)  13.0 E17-E04 S18-S11
     14 1713  911 (N07W45 L=179)   426     85        1 (1678/ 892-1744/ 939)   9.5 W41-W47 N05-N08
     15 1780  911 (N09W51 L=184)   390     44        1 (1762/ 886-1805/ 937) -56.3 W48-W53 N07-N10
     16 1828  859 (N08W58 L=191)   473     59        2 (1806/ 829-1859/ 887) -21.8 W54-W61 N06-N09
     17 1003  621 (S30W11 L=144)   233     25        1 ( 982/ 603-1026/ 642)  18.4 W09-W12 S31-S28
     18 1598  552 (S15W51 L=185)  1455    311        1 (1459/ 505-1719/ 612)  40.9 W43-W58 S21-S08
  

    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