Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1593 1528 (N27W51 L= 84)   478     93        0 (1559/1485-1624/1568) -80.5 W49-W52 N23-N29
      2  979 1416 (N31W03 L= 36)   231    126        1 ( 920/1402-1043/1441)   0.0 W00-W08 N30-N31
      3  563 1316 (N30E28 L=  5)  1632    388        5 ( 405/1267- 726/1373)   2.3 E40-E15 N27-N32
      4 1157 1185 (N14W11 L= 43)   509    134        4 (1111/1152-1226/1226) -35.5 W07-W14 N10-N16
      5 1504  837 (S13W28 L= 60)   312    122        1 (1461/ 796-1543/ 900)  45.0 W24-W30 S14-S09
      6  545  706 (S07E35 L=357)  1563    357        4 ( 498/ 572- 584/ 861)  64.8 E42-E31 S15-N02
      7  795  765 (S06E18 L= 15)   422    146        0 ( 719/ 736- 860/ 784) -26.6 E22-E13 S08-S04
      8 1640  709 (S24W37 L= 70)   686    121        2 (1587/ 668-1685/ 745)   0.0 W33-W41 S25-S21
      9  722  471 (S25E29 L=  4)   521    194        1 ( 672/ 403- 753/ 531)  45.0 E35-E26 S29-S20
     10  562  306 (S36E51 L=342)  1056    189        2 ( 520/ 252- 601/ 369)  70.0 E54-E48 S41-S30
  

    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