Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  751 1714 (N47E31 L=253)   331     42        1 ( 731/1694- 780/1733)  45.0 E33-E28 N44-N48
      2 1062 1631 (N42E02 L=282)  2870    579       10 ( 802/1513-1290/1698) -11.3 E24-W16 N35-N46
      3 1619 1508 (N36W46 L=330)  2511    354        3 (1483/1435-1776/1545)  -3.5 W32-W65 N30-N39
      4 1367 1514 (N36W23 L=307)   425     39        2 (1340/1483-1394/1541) -76.0 W20-W25 N33-N37
      5  898 1191 (N12E09 L=275)  1635    439        3 ( 865/ 982- 921/1383) -78.2 E12-E06 S01-N23
      6  267 1188 (N07E55 L=229)  2269    380        3 ( 223/1024- 304/1366) -63.4 E64-E50 S03-N17
      7 1677 1263 (N20W45 L=329)   522    129        0 (1636/1223-1702/1329)  45.0 W41-W48 N17-N24
      8 1615 1296 (N22W40 L=324)   320     66        0 (1596/1257-1634/1328) -76.0 W38-W41 N19-N24
      9  496 1289 (N15E37 L=247)   245     80        2 ( 460/1265- 542/1310)  21.8 E40-E34 N13-N16
     10  408 1209 (N09E43 L=241)   643    188        2 ( 373/1146- 447/1277) -69.4 E46-E38 N05-N13
     11  463  763 (S18E37 L=248)  1322    441        1 ( 420/ 573- 498/ 956)  72.3 E41-E33 S30-S05
     12  173  893 (S12E65 L=219)   937     87        2 ( 138/ 852- 206/ 931) -32.0 E70-E61 S14-S09
     13  894  625 (S24E06 L=278)   438    136        0 ( 838/ 589- 953/ 663)  35.5 E09-E02 S26-S20
     14 1670  537 (S26W54 L=338)  4269    463        9 (1474/ 437-1852/ 623)  21.3 W38-W77 S33-S19
     15 1320  399 (S37W28 L=312)   503    119        2 (1264/ 377-1376/ 429)  26.6 W23-W31 S39-S34
     16 1068  319 (S46W10 L=294)   378     53        1 (1036/ 302-1101/ 337)   0.0 W06-W12 S47-S44
     17  583  309 (S52E40 L=244)   735     82        1 ( 551/ 290- 625/ 335)  18.4 E45-E33 S53-S49
  

    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