Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  721 1679 (N45E31 L=188)   409     79        1 ( 684/1654- 756/1696)   0.0 E35-E28 N42-N46
      2  475 1660 (N41E53 L=166)   826    112        0 ( 431/1639- 529/1687)  29.7 E56-E48 N38-N43
      3 1620 1347 (N25W42 L=260)   512    150        0 (1574/1314-1676/1369)  24.0 W37-W46 N22-N26
      4 1804 1228 (N18W57 L=276)  1088    253        4 (1749/1147-1867/1312) -42.3 W52-W64 N12-N23
      5 1371  868 (S06W23 L=241)   458    137        1 (1343/ 799-1393/ 936) -69.4 W20-W24 S09-S01
      6 1735  588 (S23W57 L=276)   482    107        0 (1705/ 573-1776/ 601)   8.1 W53-W61 S23-S21
      7 1338  542 (S27W25 L=244)   272     95        1 (1288/ 527-1382/ 558)   0.0 W20-W27 S27-S26
      8  677  496 (S34E23 L=196)   735    134        0 ( 621/ 477- 750/ 526)  -5.7 E27-E17 S35-S32
      9  430  466 (S39E46 L=173)  1614    308        3 ( 360/ 422- 519/ 502)  12.5 E55-E37 S43-S35
  

    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