Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1519 1552 (N20W48 L=245)   801    194        2 (1445/1522-1591/1572) -18.4 W42-W54 N17-N22
      2  553 1311 (N35E23 L=174)   459    138        1 ( 494/1274- 620/1348)   5.2 E27-E16 N33-N36
      3  431 1123 (N27E35 L=162)   260     46        1 ( 412/1103- 455/1143)  14.0 E36-E32 N25-N27
      4  622  979 (N14E24 L=172)   209     50        1 ( 595/ 962- 640/1008)  26.6 E26-E22 N12-N15
      5  378  953 (N18E42 L=154)   262     62        0 ( 354/ 929- 404/ 985) -63.4 E43-E40 N16-N20
      6 1035  529 (S22E13 L=183)   222     65        2 (1007/ 503-1068/ 558) -53.1 E14-E11 S24-S19
      7  497  442 (S15E52 L=145)  1567    285        5 ( 410/ 364- 548/ 512) -76.0 E56-E46 S20-S08
  

    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