Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  823 1756 (N43E22 L= 62)   639    109        0 ( 794/1718- 868/1782)  21.8 E25-E18 N40-N45
      2  975 1703 (N40E08 L= 76)   574     75        2 ( 934/1673-1025/1727)   8.1 E11-E04 N37-N41
      3 1062 1674 (N38E01 L= 83)   507     70        1 (1028/1641-1098/1696)  45.0 E04-W01 N35-N39
      4  546 1626 (N32E39 L= 45)  1032    281        1 ( 468/1549- 599/1689) -30.5 E49-E32 N26-N36
      5 1123 1655 (N37W03 L= 87)   236     14        0 (1106/1638-1140/1674)  71.6 W01-W04 N35-N38
      6 1251 1581 (N32W13 L= 97)   819    166        1 (1174/1546-1326/1611)  -5.7 W07-W17 N29-N34
      7 1456 1520 (N29W27 L=111)   602    127        1 (1392/1497-1510/1536)  18.4 W22-W31 N27-N31
      8  318 1454 (N20E53 L= 31)   358    101        0 ( 307/1416- 330/1500) -51.3 E55-E51 N17-N22
      9 1022 1337 (N14E02 L= 82)   345     95        0 (1007/1292-1038/1394) -80.5 E02-W00 N11-N17
     10  444 1278 (N09E39 L= 45)   353     59        1 ( 424/1242- 464/1311) -53.1 E40-E37 N06-N10
     11  488  669 (S28E36 L= 48)   254     44        0 ( 471/ 642- 510/ 692)  45.0 E37-E33 S29-S26
     12 1531  548 (S30W40 L=124)  1271    415        1 (1367/ 492-1714/ 621)  21.8 W28-W53 S34-S23
  

    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