Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1591 1680 (N48W59 L=247)   564     63        0 (1560/1658-1619/1702)   7.6 W51-W66 N45-N49
      2 1206 1546 (N37W12 L=201)   298     99        1 (1166/1523-1239/1569)   9.5 W09-W15 N35-N38
      3  301 1523 (N31E67 L=121)   622     60        1 ( 275/1504- 333/1546)  -4.4 E75-E62 N30-N32
      4 1263 1488 (N33W16 L=205)   266     24        0 (1241/1467-1286/1511)   0.0 W14-W17 N31-N34
      5  485 1471 (N29E43 L=146)  1017    227        4 ( 403/1441- 607/1508)  -3.0 E50-E31 N26-N31
      6 1285 1329 (N22W17 L=205)   261     86        0 (1274/1290-1299/1375)  71.6 W15-W17 N19-N25
      7  301 1154 (N08E51 L=137)   581    138        0 ( 281/1098- 329/1221) -53.1 E54-E48 N04-N12
      8 1658 1144 (N11W43 L=232)   423     96        0 (1616/1106-1693/1175) -33.7 W39-W45 N08-N12
      9  625 1074 (N05E25 L=163)   246    109        0 ( 594/1027- 669/1102)  45.0 E27-E22 N01-N06
     10 1713  982 (N01W47 L=236)   922    181        1 (1680/ 908-1737/1056) -76.0 W44-W49 S03-N05
     11 1600  751 (S13W40 L=228)   346     66        0 (1576/ 721-1624/ 778) -31.0 W37-W42 S15-S11
     12  276  696 (S21E57 L=132)   664    130        0 ( 222/ 666- 319/ 736)  23.2 E65-E50 S23-S17
     13 1765  583 (S25W62 L=251)  2260    212        2 (1685/ 517-1818/ 665)  45.0 W54-W69 S29-S19
     14  886  580 (S25E08 L=180)   804    227        0 ( 822/ 544- 971/ 639) -16.7 E12-E02 S27-S21
     15  670  577 (S27E24 L=165)   294     79        0 ( 629/ 557- 708/ 602)  26.6 E26-E20 S28-S24
     16 1372  235 (S53W42 L=231)   737     93        1 (1336/ 214-1407/ 263) -17.4 W34-W50 S55-S50
     17 1285  216 (S56W33 L=222)   567     55        1 (1266/ 188-1301/ 245) -31.0 W28-W38 S59-S52
     18 1230  203 (S57W27 L=216)   565     60        2 (1191/ 182-1262/ 221)   0.0 W22-W31 S59-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