Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  713 1632 (N53E02 L=275)   391     57        1 ( 694/1601- 738/1665)  63.4 E04-W01 N50-N55
      2 1197 1600 (N34W31 L=307)  2057    555        6 (1019/1520-1466/1656) -43.7 W19-W44 N20-N42
      3 1618 1608 (N20W62 L=339)   444     60        0 (1584/1588-1648/1627) -36.9 W60-W64 N17-N21
      4  393 1537 (N56E41 L=235)   919    140        1 ( 350/1485- 454/1558)  -9.5 E50-E32 N51-N58
      5 1792 1453 (N06W69 L=345)   711     64        0 (1758/1424-1818/1481)   4.4 W62-W75 N04-N06
      6 1529 1414 (N13W43 L=319)   271     41        0 (1501/1395-1554/1432) -56.3 W40-W43 N11-N14
      7 1629 1311 (N04W45 L=322)   260     47        1 (1604/1283-1648/1337)  33.7 W42-W47 N03-N05
      8  823 1323 (N28E04 L=273)   245    147        0 ( 752/1300- 880/1336) -32.0 E08-W00 N25-N30
      9  280 1212 (N35E50 L=226)   703    142        0 ( 220/1183- 342/1248) -21.8 E58-E43 N32-N38
     10 1800 1063 (S15W52 L=328)   504     90        1 (1782/1017-1818/1099) -80.5 W50-W53 S17-S11
     11  915  999 (N08E07 L=270)   298    119        1 ( 876/ 967- 974/1040)  15.9 E09-E02 N06-N09
     12  344  933 (N17E46 L=230)   470    149        0 ( 299/ 909- 411/ 981) -36.9 E50-E40 N14-N20
     13  190  923 (N19E64 L=213)   858    125        0 ( 133/ 891- 240/ 949) -18.4 E72-E57 N16-N21
     14 1695  707 (S32W36 L=312)  1474    387        3 (1632/ 554-1781/ 888)  46.7 W29-W46 S41-S23
     15 1516  652 (S29W20 L=297)  1148    225        1 (1431/ 589-1573/ 751)  26.6 W13-W25 S32-S24
     16 1785  618 (S43W48 L=325)  1791    174        6 (1712/ 574-1835/ 688)  -6.3 W38-W56 S47-S38
     17  763  425 (S22E34 L=242)   281     53        0 ( 746/ 399- 785/ 450)  11.3 E36-E31 S23-S20
     18 1188  281 (S45E15 L=261)  1317    201        4 (1083/ 237-1272/ 328)  -3.2 E24-E06 S47-S43
  

    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