Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  953 1576 (N24E21 L= 22)   355     67        1 ( 921/1557- 995/1597)   9.5 E23-E17 N22-N25
      2 1558 1519 (N41W21 L= 64)   392     50        1 (1527/1501-1587/1546)  21.8 W18-W23 N39-N42
      3  407 1436 (N02E51 L=352)   349     29        0 ( 383/1415- 433/1456)  33.7 E53-E49 N00-N03
      4 1808 1125 (N24W50 L= 93)  2275    520        5 (1688/ 932-1894/1344) -33.7 W35-W63 N14-N34
      5  592  980 (S19E24 L= 19)   400    166        0 ( 521/ 950- 668/1005)  33.7 E28-E19 S22-S16
      6  434  875 (S30E34 L=  9)   298     84        0 ( 401/ 846- 472/ 895)  31.0 E36-E31 S31-S27
      7  237  814 (S38E55 L=348)   532     87        1 ( 191/ 801- 282/ 837)   8.1 E63-E49 S39-S36
      8 1488  298 (S31W65 L=108)   871    103        2 (1457/ 264-1528/ 345)  -7.1 W58-W74 S31-S28
      9  961  295 (S51W27 L= 70)   284     47        0 ( 942/ 271- 982/ 320)  -8.1 W23-W30 S52-S49
  

    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