Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  910 1684 (N34E28 L=265)   978    136        1 ( 855/1638- 965/1725) -18.4 E33-E21 N32-N36
      2  568 1651 (N24E51 L=242)   439     64        0 ( 551/1615- 586/1683) -53.1 E53-E48 N21-N25
      3 1506 1583 (N46W21 L=314)   338     78        0 (1482/1556-1541/1601)  18.4 W18-W24 N42-N47
      4  982 1382 (N17E11 L=282)   215     49        1 ( 960/1362-1009/1402)  45.0 E12-E09 N14-N17
      5  223 1206 (S10E62 L=231)  2563    195        8 ( 182/1120- 282/1310) -42.3 E71-E55 S14-S04
      6  144 1134 (S16E72 L=221)  1633    118        3 ( 118/1084- 183/1185)  37.9 E78-E66 S18-S11
      7 1893 1139 (N27W65 L=358)   948     84        2 (1861/1091-1911/1184)  60.3 W59-W70 N23-N30
      8 1515 1085 (N12W28 L=321)   285    108        0 (1486/1030-1540/1126) -45.0 W25-W30 N09-N14
      9 1656  546 (S14W58 L=350)  1471    251        6 (1601/ 471-1698/ 619) -36.0 W52-W64 S18-S10
     10 1428  494 (S23W42 L=335)   366     75        0 (1416/ 458-1443/ 545) -51.3 W39-W44 S25-S20
     11 1544  473 (S21W52 L=345)  1291    235        3 (1447/ 411-1626/ 525)  45.0 W46-W56 S26-S16
     12 1054  220 (S54W37 L=330)   488     53        0 (1026/ 196-1079/ 246) -11.3 W32-W42 S55-S52
  

    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