Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  701 1675 (N24E42 L=101)   376     72        0 ( 676/1637- 716/1707) -29.7 E46-E39 N21-N25
      2  404 1644 (N15E67 L= 76)  2433    255        9 ( 331/1563- 462/1693)  24.8 E76-E57 N10-N19
      3  538 1615 (N16E51 L= 92)   343     53        2 ( 518/1592- 558/1643) -63.4 E52-E48 N14-N18
      4  382 1601 (N12E65 L= 78)   488     52        0 ( 366/1573- 403/1622)  90.0 E68-E62 N10-N13
      5  675 1579 (N17E39 L=105)   307     83        1 ( 635/1555- 702/1603)  26.6 E41-E36 N16-N19
      6 1718 1451 (N43W39 L=182)   426     46        0 (1692/1433-1739/1478)  45.0 W35-W42 N41-N46
      7 1390  547 (S21W38 L=181)   613    205        0 (1328/ 486-1449/ 608)  84.3 W35-W39 S25-S15
      8  956  534 (S35W12 L=155)   219     56        0 ( 925/ 521- 989/ 547)  26.6 W09-W13 S36-S33
      9  562  487 (S52E18 L=125)   555     70        0 ( 522/ 461- 607/ 511)  -5.7 E22-E12 S53-S50
  

    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