Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  854 1736 (N32E32 L=314)   543    114        1 ( 801/1695- 884/1780) -21.8 E37-E27 N29-N34
      2 1007 1573 (N24E15 L=331)   291     75        1 ( 987/1539-1026/1616) -51.3 E16-E12 N22-N27
      3  880 1530 (N19E21 L=325)   442    158        0 ( 805/1507- 961/1551)  11.3 E25-E15 N17-N20
      4 1606 1358 (N27W29 L= 15)   476     96        1 (1574/1312-1649/1388)  54.5 W26-W31 N23-N30
      5 1705 1291 (N26W37 L= 23)   335     92        0 (1673/1247-1736/1326) -26.6 W34-W40 N23-N27
      6 1322 1280 (N14W11 L=357)   343    116        2 (1300/1232-1365/1313) -78.7 W09-W13 N11-N16
      7  976 1222 (N03E07 L=339)   615    259        0 ( 850/1196-1089/1257)  24.8 E14-E01 N00-N07
      8 1608 1204 (N17W31 L= 17)   288     58        2 (1584/1181-1633/1226)  53.1 W29-W32 N15-N19
      9 1857  988 (N13W57 L= 43)   840    147        0 (1844/ 922-1869/1067) -76.0 W54-W58 N09-N17
     10  213  971 (S26E57 L=289)  1363    192        4 ( 121/ 938- 291/1011)   9.9 E71-E48 S28-S23
     11 1528  919 (S01W32 L= 18)   207     51        0 (1505/ 900-1553/ 937)   0.0 W29-W33 S01-N00
     12 1478  719 (S12W35 L= 21)   266     40        2 (1449/ 696-1508/ 742)  33.7 W33-W36 S14-S10
     13  202  668 (S45E67 L=279)  1244    121        2 ( 144/ 640- 261/ 692)  -7.4 E83-E56 S46-S43
     14 1173  638 (S25W19 L=  5)   305     59        0 (1137/ 615-1204/ 659)   0.0 W16-W21 S25-S23
     15 1344  640 (S20W29 L= 15)   219     98        0 (1292/ 630-1392/ 648)  15.9 W25-W32 S21-S18
     16 1464  500 (S25W44 L= 30)  2060    483        3 (1261/ 447-1703/ 568)  36.5 W32-W59 S33-S13
  

    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