Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1565 1644 (N33W48 L=292)   378     42        0 (1539/1630-1588/1660) -14.0 W45-W51 N31-N33
      2  309 1415 (N30E49 L=195)  1085    281        2 ( 247/1375- 383/1465) -28.1 E56-E41 N26-N34
      3 1185 1437 (N22W14 L=258)   558    216        0 (1076/1416-1294/1458)  -4.1 W07-W21 N21-N23
      4 1751 1276 (N07W52 L=296)   917    222        1 (1677/1240-1807/1333) -15.3 W46-W58 N05-N11
      5  830 1310 (N18E09 L=234)   241     60        1 ( 796/1296- 867/1327)   0.0 E11-E06 N16-N18
      6  359 1248 (N19E43 L=201)   282     56        0 ( 341/1222- 379/1278)  56.3 E44-E40 N17-N20
      7  160 1108 (N12E62 L=181)  2811    402        7 ( 105/1005- 220/1214)  42.5 E71-E55 N06-N19
      8 1346 1070 (S01W20 L=263)   528    156        0 (1285/1035-1417/1123) -35.5 W16-W23 S03-N02
      9  173  985 (N05E62 L=182)   488     78        1 ( 152/ 953- 199/1013) -45.0 E64-E59 N02-N06
     10 1181  862 (S12W08 L=252)   326     56        1 (1153/ 830-1204/ 892)  45.0 W06-W09 S13-S10
     11 1822  789 (S22W56 L=299)  1423    121        4 (1784/ 719-1855/ 873)  66.0 W51-W58 S26-S17
     12 1181  689 (S23W07 L=251)   519    198        0 (1118/ 595-1260/ 780)  45.0 W01-W12 S28-S17
     13 1051  696 (S21E01 L=242)   231     94        0 ( 992/ 667-1078/ 730)  31.0 E05-W00 S22-S19
     14  174  628 (S16E72 L=171)   785     93        2 ( 158/ 585- 193/ 663) -68.2 E75-E69 S18-S13
     15  338  607 (S19E54 L=190)   610    153        2 ( 309/ 572- 366/ 639) -53.1 E56-E49 S21-S17
     16  360  507 (S26E55 L=188)   848    166        0 ( 326/ 467- 406/ 547)  16.7 E59-E49 S28-S23
     17 1324  376 (S47W17 L=261)   373     56        1 (1297/ 351-1347/ 404)  53.1 W14-W19 S48-S44
     18  966  340 (S45E14 L=230)   505     93        1 ( 933/ 310- 991/ 386)  39.8 E16-E10 S47-S41
  

    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