Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1424 1789 (N45W63 L=352)  2399    181        5 (1373/1724-1462/1855)  10.0 W54-W88 N39-N48
      2  792 1742 (N60E03 L=286)   602     82        1 ( 771/1702- 819/1783)  49.4 E05-W01 N56-N63
      3 1332 1552 (N32W34 L=323)  1982    353        5 (1218/1477-1517/1624) -45.0 W25-W44 N23-N39
      4  421 1439 (N42E41 L=248)  1823    360        4 ( 267/1395- 566/1489) -10.4 E64-E26 N37-N46
      5 1833 1378 (N09W70 L=359)   565     70        1 (1821/1336-1846/1412)  33.7 W65-W74 N07-N11
      6 1053 1084 (N10W03 L=292)   532    183        2 ( 992/1017-1125/1145) -60.9 W00-W05 N04-N14
      7 1848 1029 (S10W59 L=348)   571    105        0 (1831/ 983-1869/1077)  33.7 W56-W62 S12-S07
      8 1320  813 (S11W14 L=303)   612    285        7 (1275/ 717-1371/ 902) -90.0 W11-W17 S17-S04
      9 1435  858 (S11W22 L=311)   236     94        0 (1410/ 812-1468/ 892)  45.0 W20-W24 S13-S08
     10 1586  608 (S30W31 L=320)  1987    432        3 (1434/ 513-1721/ 768)  15.6 W18-W43 S35-S22
     11  445  473 (S19E53 L=237)  1282    274        7 ( 376/ 376- 529/ 578) -82.4 E55-E48 S26-S11
     12  885  232 (S44E32 L=257)  2833    431        4 ( 732/ 171-1068/ 306)   7.1 E52-E12 S48-S41
  

    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