Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  834 1821 (N62W19 L=117)   526    104        1 ( 802/1787- 882/1864)  10.6 W13-W30 N59-N64
      2 1144 1787 (N46W39 L=137)   502     83        1 (1117/1745-1173/1814)  35.5 W36-W43 N43-N48
      3 1373 1730 (N35W51 L=148)  1096    134        2 (1307/1702-1425/1758) -35.5 W44-W54 N31-N37
      4 1465 1640 (N27W50 L=147)   560     88        1 (1444/1601-1491/1681)  51.3 W47-W52 N24-N29
      5  774 1595 (N45W00 L= 98)   254     62        1 ( 749/1563- 797/1621) -90.0 E01-W01 N42-N47
      6  321 1484 (N51E49 L= 48)   394     65        0 ( 301/1454- 342/1513) -26.6 E55-E45 N48-N53
      7  955 1357 (N25W05 L=102)   328    107        1 ( 931/1316- 985/1407) -90.0 W02-W06 N21-N27
      8  530 1385 (N38E24 L= 73)   320     89        1 ( 493/1355- 569/1404) -45.0 E27-E21 N34-N39
      9  247 1381 (N45E57 L= 40)   706     75        1 ( 214/1351- 276/1399) -18.4 E64-E51 N42-N46
     10  220 1060 (N25E56 L= 41)  3424    577        8 ( 172/ 842- 284/1293) -74.1 E64-E50 N11-N39
     11  479 1225 (N29E30 L= 67)   234     84        1 ( 454/1203- 514/1247) -11.3 E32-E27 N27-N30
     12  345 1108 (N25E43 L= 55)   280     71        1 ( 323/1079- 381/1146)  63.4 E45-E39 N23-N27
     13 1254 1089 (N02W14 L=112)   290     63        1 (1226/1060-1279/1121)  14.0 W12-W16 N00-N03
     14 1129  828 (S09W01 L= 98)   514    165        0 (1105/ 747-1163/ 893)  49.4 E02-W04 S13-S05
     15 1750  645 (S39W44 L=141)   759    102        2 (1704/ 619-1807/ 665)  -4.4 W38-W51 S40-S36
     16  902  464 (S25E23 L= 74)   396    131        0 ( 833/ 447- 960/ 479) -14.0 E27-E19 S26-S23
     17  783  382 (S27E35 L= 63)   273     96        0 ( 764/ 347- 819/ 416)  33.7 E36-E30 S29-S25
     18 1229  241 (S53E13 L= 84)   946    125        3 (1166/ 217-1289/ 272) -26.6 E17-E07 S55-S49
     19 1125  166 (S56E29 L= 68)   468     56        0 (1100/ 149-1147/ 194)  12.5 E32-E23 S58-S53
  

    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