Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1502 1717 (N31W59 L=232)  1070    202        0 (1393/1696-1578/1740) -20.6 W49-W66 N28-N34
      2 1122 1662 (N39W26 L=199)   980    229        0 (1008/1618-1217/1700)  -3.2 W16-W34 N37-N39
      3 1323 1447 (N20W29 L=202)   259     74        1 (1285/1425-1359/1469) -53.1 W27-W30 N18-N22
      4  358 1286 (N33E41 L=132)   944    189        1 ( 267/1264- 444/1309) -14.7 E51-E32 N30-N36
      5  867 1192 (N16E05 L=168)   183     53        0 ( 834/1180- 901/1204) -33.7 E07-E03 N15-N17
      6  536 1164 (N22E27 L=146)   309     84        0 ( 489/1146- 579/1179)   0.0 E30-E23 N20-N22
      7 1668 1030 (S12W40 L=213)   477    176        1 (1638/ 945-1702/1109)  45.0 W36-W44 S15-S07
      8  181  976 (N18E62 L=111)  2151    250        4 (  84/ 918- 259/1052)   5.5 E82-E51 N14-N21
      9 1781 1006 (S16W49 L=223)   473     74        0 (1766/ 965-1794/1044)  68.2 W47-W51 S18-S13
     10 1727  989 (S15W44 L=217)   290     61        0 (1709/ 956-1748/1025)  45.0 W42-W46 S17-S13
     11 1118  952 (S03W04 L=177)   289    110        0 (1060/ 939-1171/ 980)  -8.1 W00-W07 S03-S01
     12 1839  901 (S24W55 L=228)  1172    133        2 (1799/ 855-1869/ 972)  50.2 W50-W58 S26-S19
     13 1400  733 (S22W16 L=189)   267     55        0 (1387/ 694-1414/ 768)  53.1 W14-W17 S24-S20
     14 1204  626 (S24W01 L=174)   807    252        2 (1149/ 520-1261/ 741)  47.7 E04-W07 S29-S18
     15 1009  606 (S20E11 L=162)   864    254        3 ( 966/ 485-1058/ 723)  49.8 E16-E05 S27-S14
     16  691  479 (S21E35 L=138)  1426    395        1 ( 639/ 305- 752/ 645)  39.6 E48-E25 S31-S11
     17 1263  536 (S31W03 L=176)   375    120        0 (1205/ 508-1311/ 574)   8.1 E01-W06 S32-S29
     18 1513  400 (S50W21 L=194)  1090    197        0 (1411/ 362-1588/ 439)  -3.0 W10-W30 S52-S47
  

    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