Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  404 1423 (N29E42 L= 17)  4177    502        7 ( 266/1233- 520/1618) -30.8 E78-E30 N16-N45
      2 1572 1488 (N22W43 L=102)   766    174        1 (1483/1459-1652/1513) -17.1 W36-W49 N19-N24
      3  163 1182 (N17E64 L=354)  4578    397        8 (  86/1015- 251/1294)  41.4 E78-E54 N07-N23
      4 1489  966 (S09W28 L= 86)   274     53        1 (1465/ 940-1507/ 998)  76.0 W26-W29 S10-S06
      5  449  706 (S15E42 L= 17)  1049    266        2 ( 372/ 600- 523/ 805)  32.9 E51-E34 S20-S09
      6  250  617 (S18E65 L=353)  2540    358        2 ( 212/ 468- 302/ 791)  63.4 E72-E58 S27-S07
      7 1293  744 (S20W14 L= 73)   275    115        1 (1255/ 714-1353/ 767)  15.9 W11-W18 S21-S19
      8 1108  672 (S23W02 L= 60)   810    266        4 (1004/ 638-1205/ 744) -36.0 E03-W08 S26-S17
      9 1485  464 (S41W31 L= 89)   842    101        2 (1438/ 435-1530/ 506)   0.0 W27-W35 S43-S38
  

    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