Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  768 1828 (N64E36 L=322)   451     52        0 ( 741/1809- 797/1844)  31.0 E39-E31 N61-N65
      2  960 1823 (N63E07 L=351)   428     68        1 ( 933/1795- 979/1843) -29.7 E10-E03 N60-N65
      3 1659 1502 (N32W55 L= 52)   385    124        0 (1618/1454-1698/1548) -66.8 W52-W57 N28-N35
      4  351 1440 (N30E56 L=302)   580    114        1 ( 307/1411- 387/1467)  20.6 E59-E50 N27-N31
      5 1701 1386 (N24W54 L= 51)  1203    246        4 (1585/1314-1780/1459) -25.2 W44-W61 N19-N29
      6 1358 1374 (N25W24 L= 21)   376    177        0 (1283/1318-1415/1405) -33.7 W18-W27 N21-N27
      7  309 1330 (N22E56 L=302)   560    117        1 ( 284/1294- 352/1382) -26.6 E60-E49 N20-N25
      8 1559 1340 (N22W39 L= 36)   300     76        1 (1521/1315-1583/1369)  21.8 W35-W40 N20-N24
      9  821 1307 (N22E13 L=344)   255     82        0 ( 794/1273- 838/1346)  68.2 E14-E12 N19-N24
     10 1843  805 (S13W64 L= 62)  1087    184        3 (1819/ 732-1870/ 894)  84.3 W60-W68 S17-S07
     11 1748  859 (S09W52 L= 49)   535     46        1 (1718/ 829-1772/ 884)  56.3 W48-W54 S10-S07
     12 1778  774 (S14W56 L= 54)   810     80        4 (1745/ 743-1814/ 798)   7.1 W52-W60 S16-S12
     13 1458  699 (S18W29 L= 26)   637    172        1 (1394/ 655-1536/ 738)  29.1 W24-W34 S20-S15
     14  512  673 (S18E36 L=321)   354     98        0 ( 495/ 633- 526/ 727)  66.8 E37-E34 S21-S14
     15 1797  612 (S26W66 L= 64)  1140    120        1 (1741/ 586-1848/ 643)   0.0 W57-W77 S27-S24
     16  878  415 (S37E12 L=345)  2847    324        8 ( 791/ 333-1004/ 537) -13.2 E19-E02 S43-S27
     17 1086  385 (S40W04 L=  1)   451     55        2 (1058/ 359-1115/ 415) -45.0 W01-W06 S41-S37
     18  601  310 (S47E43 L=314)   626     84        0 ( 589/ 266- 614/ 354)  39.3 E49-E38 S51-S42
  

    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