Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1556 1687 (N25W59 L=182)  2017    251        3 (1460/1639-1631/1716) -39.3 W52-W65 N20-N30
      2  751 1332 (N30E08 L=115)   855    369        2 ( 583/1300- 943/1369) -19.2 E19-W04 N25-N35
      3 1696 1267 (S01W48 L=171)   659    109        3 (1659/1239-1735/1303) -21.8 W44-W51 S02-N01
      4 1759 1237 (S04W53 L=175)   313     54        1 (1741/1213-1781/1261)  45.0 W50-W54 S05-S02
      5 1635 1004 (S13W36 L=158)  1379    431        2 (1600/ 821-1698/1223)  64.4 W31-W46 S24-S01
      6 1020 1189 (N14W04 L=127)   390     65        0 ( 988/1157-1052/1222) -63.4 W03-W05 N12-N16
      7 1848 1080 (S16W56 L=178)   622    121        1 (1806/1050-1877/1114)  -7.1 W51-W59 S17-S14
      8 1448  994 (S08W23 L=146)   271     39        2 (1422/ 970-1469/1024)  33.7 W21-W24 S09-S06
      9 1328  788 (S16W11 L=134)   523    144        1 (1251/ 749-1384/ 821) -53.1 W07-W13 S19-S11
     10  855  504 (S20E25 L= 98)   325     81        1 ( 818/ 481- 878/ 539)  21.8 E27-E22 S21-S17
     11  502  359 (S21E56 L= 67)   476     86        2 ( 484/ 328- 520/ 387)  38.7 E58-E53 S22-S18
     12 1196  275 (S47E15 L=108)  1361    188        3 (1150/ 226-1239/ 331)  63.4 E20-E10 S51-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