Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  787 1761 (N53E28 L=160)   324     33        1 ( 765/1742- 804/1779)  90.0 E30-E25 N51-N55
      2 1594 1681 (N48W60 L=248)   549     47        0 (1563/1658-1621/1701)   8.1 W52-W66 N45-N49
      3 1212 1547 (N37W13 L=201)   290     64        1 (1181/1523-1242/1568)  31.0 W10-W15 N35-N38
      4  302 1523 (N31E67 L=121)   619     63        1 ( 276/1504- 335/1548)  -4.4 E75-E62 N30-N32
      5 1266 1488 (N33W16 L=205)   265     39        0 (1243/1465-1289/1509)   0.0 W14-W18 N31-N34
      6  488 1470 (N29E42 L=146)  1028    230        3 ( 410/1440- 607/1508)  -3.0 E50-E31 N26-N31
      7 1288 1326 (N22W17 L=205)   276    100        0 (1272/1274-1304/1370)  74.1 W15-W17 N18-N25
      8  302 1154 (N08E51 L=137)   614    132        1 ( 281/1096- 331/1223) -53.1 E54-E47 N04-N12
      9 1660 1145 (N11W43 L=232)   464    115        0 (1612/1103-1696/1177) -33.7 W39-W45 N08-N13
     10  625 1072 (N05E25 L=163)   234    117        0 ( 594/1023- 670/1100)  45.0 E27-E22 N01-N06
     11 1716  979 (N01W47 L=236)   895    174        1 (1683/ 909-1739/1058)  90.0 W44-W49 S03-N05
     12 1601  752 (S13W40 L=228)   343     53        0 (1578/ 724-1622/ 777) -56.3 W37-W41 S14-S11
     13  279  697 (S21E56 L=132)   610    126        0 ( 230/ 666- 321/ 737)  21.0 E63-E50 S22-S17
     14 1765  583 (S25W62 L=251)  2079    209        2 (1686/ 518-1820/ 665)  48.4 W54-W68 S29-S19
     15  888  581 (S25E08 L=181)   847    230        2 ( 824/ 544- 978/ 639) -15.3 E12-E01 S27-S21
     16  673  576 (S27E23 L=165)   295     80        0 ( 629/ 559- 709/ 601)  26.6 E26-E20 S27-S24
     17 1373  235 (S53W43 L=231)   633     93        2 (1339/ 217-1408/ 262) -18.4 W35-W50 S55-S50
     18 1287  215 (S56W34 L=222)   604     52        2 (1268/ 189-1308/ 245) -33.7 W28-W38 S58-S52
     19 1218  202 (S57W26 L=215)   658     89        2 (1162/ 182-1258/ 220)   5.2 W19-W30 S59-S55
  

    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