Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  527 1545 (N41E39 L= 64)   339     54        0 ( 513/1514- 541/1577) -68.2 E41-E36 N38-N43
      2 1029 1557 (N39W04 L=107)   381     90        1 ( 979/1533-1072/1569)   0.0 W00-W07 N37-N40
      3 1238 1423 (N28W18 L=121)   351    162        0 (1207/1350-1258/1487) -90.0 W15-W19 N23-N32
      4  509 1308 (N25E35 L= 68)   478     98        0 ( 495/1258- 526/1354) -63.4 E36-E33 N21-N27
      5 1818 1245 (N11W62 L=165)   588     78        2 (1799/1198-1841/1281) -71.6 W60-W64 N07-N13
      6 1778 1242 (N11W57 L=160)   367     73        0 (1756/1206-1797/1278) -90.0 W55-W58 N09-N13
      7  196  881 (S02E64 L= 39)  1323    243        2 ( 177/ 773- 225/ 988)  63.4 E67-E58 S08-N04
      8 1778  765 (S18W55 L=157)  4595    565       19 (1717/ 564-1857/ 924) -82.9 W47-W64 S31-S07
      9 1396  803 (S11W22 L=125)   258     90        0 (1357/ 770-1429/ 839)  45.0 W19-W24 S13-S09
     10  402  655 (S16E47 L= 56)   294     88        0 ( 371/ 615- 419/ 692) -68.2 E48-E45 S18-S13
     11  863  463 (S31E15 L= 87)  1691    440        5 ( 685/ 393- 989/ 603)  23.7 E30-E05 S35-S21
     12  420  495 (S27E51 L= 51)   386     88        2 ( 392/ 479- 474/ 511)   0.0 E55-E45 S28-S25
  

    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