Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  829 1817 (N62W19 L=117)   434     99        0 ( 797/1784- 872/1853)  11.3 W12-W27 N60-N64
      2 1157 1783 (N46W40 L=138)   756    159        2 (1108/1744-1221/1805) -10.3 W34-W46 N43-N48
      3 1363 1724 (N35W50 L=148)   926    139        3 (1301/1694-1415/1750) -50.2 W44-W53 N31-N37
      4  925 1683 (N47W16 L=114)   282     60        2 ( 885/1665- 948/1705) -21.8 W12-W18 N45-N48
      5 1454 1632 (N27W49 L=147)   409     74        1 (1433/1595-1477/1669)  59.0 W45-W51 N24-N29
      6  940 1353 (N25W04 L=102)   304    112        2 ( 918/1309- 972/1404) -90.0 W01-W05 N21-N28
      7  520 1384 (N38E25 L= 73)   253     84        0 ( 485/1359- 556/1400) -38.7 E28-E22 N35-N39
      8  217 1143 (N30E57 L= 41)  3924    716       18 ( 157/ 888- 279/1398) -86.4 E66-E50 N14-N46
      9  795 1322 (N27E06 L= 92)   261    107        1 ( 759/1276- 819/1362) -71.6 E07-E05 N23-N29
     10  484 1244 (N30E30 L= 69)   292     96        1 ( 454/1202- 506/1288)  45.0 E32-E27 N27-N32
     11 1102 1247 (N15W10 L=108)   220    126        2 (1065/1220-1144/1272) -11.3 W07-W12 N13-N16
     12  232 1105 (N27E55 L= 43)   367     67        2 ( 200/1084- 256/1148)  45.0 E59-E52 N25-N29
     13 1236 1072 (N02W13 L=111)   196     32        0 (1215/1052-1254/1091) -71.6 W11-W14 N00-N03
     14 1116  819 (S09W00 L= 98)   498    191        0 (1089/ 720-1149/ 893)  56.3 E03-W03 S14-S05
     15 1744  642 (S39W43 L=142)   708    101        2 (1699/ 618-1797/ 662)  -4.8 W37-W50 S41-S37
     16  900  464 (S25E23 L= 75)   371    129        0 ( 832/ 450- 958/ 478) -20.6 E27-E19 S26-S23
     17 1241  242 (S53E12 L= 86)   539     53        1 (1217/ 216-1270/ 272)  74.1 E14-E08 S56-S49
     18 1121  166 (S57E30 L= 69)   446     48        0 (1095/ 148-1144/ 187) -11.3 E33-E25 S58-S54
  

    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