Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  565 1646 (N60E18 L= 77)   323     39        0 ( 544/1624- 591/1668)   9.5 E21-E14 N58-N61
      2  748 1608 (N50W00 L= 95)   318     64        1 ( 724/1580- 767/1640)  63.4 E02-W02 N47-N51
      3 1583 1586 (N20W57 L=152)   800     86        3 (1541/1548-1618/1623) -90.0 W54-W60 N17-N22
      4 1660 1587 (N17W65 L=160)   668    102        3 (1619/1547-1690/1622) -68.2 W61-W72 N14-N19
      5  352 1327 (N42E43 L= 52)   593    124        1 ( 336/1264- 365/1395) -77.5 E45-E40 N37-N46
      6  661 1142 (N22E19 L= 76)   784    279        1 ( 548/1045- 762/1198) -51.8 E25-E13 N14-N28
      7  338 1103 (N28E45 L= 50)   603    137        0 ( 295/1052- 372/1166)  32.0 E49-E40 N25-N30
      8 1765 1065 (S13W48 L=143)  1879    333        3 (1616/ 984-1863/1152)   0.0 W34-W61 S15-S10
      9  366  993 (N20E44 L= 51)  1051    163        3 ( 323/ 937- 405/1048)  26.6 E47-E39 N16-N23
     10  191  948 (N21E63 L= 32)   542     62        0 ( 164/ 914- 219/ 981) -33.7 E66-E59 N18-N23
     11 1522  693 (S27W21 L=116)  1187    309        1 (1454/ 581-1575/ 840)  42.7 W14-W27 S33-S19
     12 1276  598 (S25W03 L= 98)  1658    578        0 (1113/ 420-1468/ 791)  22.8 E12-W19 S32-S18
     13 1693  718 (S31W36 L=131)   336     39        1 (1676/ 693-1713/ 744)  56.3 W34-W38 S32-S29
     14 1826  666 (S41W54 L=149)   813     92        0 (1782/ 640-1872/ 688)  -7.1 W46-W64 S43-S40
     15  604  379 (S21E48 L= 47)   615     70        2 ( 571/ 348- 639/ 405)  12.5 E53-E44 S22-S19
     16 1518  340 (S56W17 L=111)   415     37        1 (1495/ 322-1542/ 364)  11.3 W13-W19 S57-S53
  

    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